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APRESENTAÇÃO 


O e-book: “Engenharia sanitária e ambiental: Recursos hídricos e 
tratamento de água 3” é constituído por cinco capítulos de livro que tratam da 
disponibilidade, qualidade e principais uso de recursos hídricos para fins potáveis 
ou não por todos os segmentos da sociedade. 

O primeiro capítulo apresenta um estudo no qual se discute a importância 
do monitoramento constante em relação à segurança de barragens para 
armazenamento de rejeitos provenientes de atividades de mineração e/ou 
recursos hídricos, bem como a descrição de todas as legislações nacionais 
existentes e vigentes no território brasileiro. O capítulo 2 se propôs a apresentar 
um estudo de comparação de cálculos de vazões por meio do uso do Perfilador 
Acústico de Corrente por efeito Doppler (ADCP) em modo estático e o método de 
máxima entropia M em rios brasileiros monitorados pela Rede Hidrometeorológica 
Nacional — RHN presentes no estado da Bahia e Sergipe. 

O terceiro capítulo avaliou a importância da captação de água de chuva, 
bem como a sua utilização para reduzir à necessidade básica de populações que 
não possuem acesso a água para fins potáveis ou não, bem como a redução 
de etapas de tratamento de água que geraria economia para a população 
beneficiada. O capítulo 4 analisou os diferentes impactos provenientes do 
Projeto de Integração do rio São Francisco (PISF) em relação ao canal Acauã- 
Araçagi, bem como estimar as tarifas de água a ser cobrada pela população 
beneficiária deste adutor com vazão máxima de 10 mº/s e uma extensão de 112 
km, utilizando-se testes simulatórios com o ModSIM P32. 

Por fim, o quinto capítulo apresenta uma proposta de utilização da argila 
como adsorvente para a forma mais tóxica do arsênio presente em concentrações 
traços e/ou ultra-traços em diferentes compartimentos aquáticos a partir da 
aplicação de um dispositivo de extração miniaturizado em ponteira descartável. 

Nesta perspectiva, a Atena Editora vem trabalhando de forma a estimular 
e incentivar cada vez mais pesquisadores do Brasil e de outros países a 
publicarem seus trabalhos com garantia de qualidade e excelência em forma de 
livros, capítulos de livros e artigos científicos. 


Cleiseano Emanuel da Silva Paniagua 
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CAPÍTULO 1 


ÍNDICE SUSTENTÁVEL AMBIENTAL DE 
SEGURANÇA DE BARRAGENS - ISASB 


Maria Bernardete Guimarães 

Analista de Meio Ambiente e Recursos 
hídricos do IEMA-ES-Brasil; Engenheira 
Ambiental,M.Sc. Eng.º? Ambiental 


PALAVRAS-CHAVE: Índice 


Barragens; 


Segurança 
Acidentes com Barragens; 
Índice Sustentável. 


INTRODUÇÃO 


A lei 12.334 de 2010 estabeleceu 
que a Política Nacional de Segurança 
de Barragens, alterada pela Lei nº 
14.066 de 2020, aplica-se a barragens 
destinadas a: acumulação de água para 
quaisquer usos; à disposição temporária 
de rejeitos; à acumulação de resíduos 
industriais, que apresentem pelo menos 
uma das seguintes características: altura 
do maciço, medida do encontro do pé do 
talude de jusante com nível do solo até a 
crista de coroamento do barramento, maior 
ou igual a 15 (quinze) metros (redação da 
lei 14.066 de 2020); capacidade total do 
reservatório maior ou igual a 3.000.000 
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metros cúbicos; reservatório que contenha 
resíduos 
técnicas aplicáveis (ABNT NBR 10004 de 
2004-Classificação de Resíduos Sólidos); 


perigosos conforme normas 


categoria de dano potencial associado 
médio ou alto, em termos econômicos, 
sociais, ambientais ou de perda de vidas 
humanas, conforme definido no art. 7º 
desta Lei (redação dada pela Lei nº 
14.066 de 2020); categoria de risco alto, 
a critério do órgão fiscalizador, conforme 
definido no art. 7º desta Lei (incluído na Lei 
14.066 de 2020). NEVES (2018) mostra as 
legislações comentadas sobre Segurança 
de Barragens. 

Barragem é qualquer estrutura 
construída dentro ou fora de um curso 
de água permanente ou temporário, em 
talvegue ou em cava exaurida com dique, 
para fins de contenção ou acumulação 
de susbtâncias líquidas ou de misturas 
de líquidos e sólidos, compreendendo o 
barramento e as estruturas associadas 
(BRASIL, 2020). São instrumentos da 
Política 


Nacional de Segurança de 


Barragens- PNSBo sistema de classificação 
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de barragens por categoria de risco e por dano potencial associado; o plano de segurança 
de barragens; o sistema nacional de informações sobre segurança de barragens- SNISB; o 
sistema nacional de informação sobre Meio Ambiente- SINIMA; o cadastro técnico federal 
de atividades e instrumentos de defesa ambiental; o cadastro técnico federal de atividades 
potencialmente poluidoras ou utilizadoras de recursos ambientais; o relatório de segurança 
de barragens (BRASIL, 2010). O dano potencial associado à barragem é o dano que pode 
ocorrer devido a rompimentos, vazamentos, infiltração no solo ou mau funcionamento de 
uma barragem, independente da sua probabilidade de ocorrência a ser graduado de acordo 
com as perdas de vidas humanas e os impactos sociais, econômicos e ambientais (BRASIL, 
2020). O Plano de Segurança de Barragens-PSB é um instrumento da PNSB- Política 
Nacional de Segurança de Barragens, cujo objetivo principal é auxiliar o empreendedor na 
garantia de padrões de segurança da Barragem reduzindo os riscos de acidentes e suas 
consequências, incluindo, quando aplicável, o Plano de Ação de Emergência-PAE. 

O Plano também auxilia as entidades fiscalizadoras de segurança de barragem 
na verificação do cumprimento das normas concernentes ao tema. Deve conter dados 
técnicos da barragem, de construção, operação, manutenção e panorama do estado atual 
da segurança por meio das inspeções realizadas, devendo principalmente, servir como 
uma ferramenta de planejamento de gestão da segurança da Barragem. 

A nova legislação, Lei nº 14.006 de 2020 traz como alterações à lei 12,334 de 2010 
os seguintes itens: A PNSB deverá estabelecer programa de educação ambiental e de 
comunicação sobre segurança de barragem, com o objetivo de conscientizar a sociedade 
da importância da segurança de barragens e de desenvolver a cultura de prevenção 
a acidentes e desastres, que deverá contemplar as seguintes medidas: fica proibida a 
construção ou o alteamento de barragem de mineração pelo método a montante (entende- 
se por alteamento a montante a metodologia construtiva de barragem em que os diques 
de contenção se apoiam sobre seu próprio rejeito ou sedimento previamente lançado 
depositado; o empreendedor deve concluir a descaracterização da barragem construída 
ou alterada pelo método a montante até 25 de fevereiro de 2022, considerada a solução 
técnica exigida pela entidade que regula e fiscaliza a atividade minerária e pela autoridade 
licenciadora do sistema nacional do meio ambiente (SISNAMA). O artigo 18 veta a 
implantação de barragem de mineração cujos estudos de cenários de ruptura identificaram 
a existência de comunidade na ZAS. Na figura 1 a Gestão de Riscos, de acordo com a ISO 
31.000 (ABNT,2011). 
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Figura 1- ABNT (2011)- ISO 31.000 Gestão de Riscos 
FONTE: ABNT (2011) EDITADO PELA AUTORA 


Quanto àfiscalização a Lei 14.006 de 2020 diz que os órgãos fiscalizadores devemcriar 
sistema de credenciamento de pessoas físicas e jurídicas habilitadas a atestar a segurança 
da barragem, incluída a certificação, na forma do regulamento. O laudo técnico referente 
às causas do rompimento de barragem deve ser elaborado por peritos independentes, às 
expensas do empreendedor, em coordenação do órgão fiscalizador. Também diz que o 
poder público municipal deve adotar as medidas necessárias para impedir o parcelamento, 
o uso e a ocupação do solo urbano na ZAS- Zona de Altosalvamento. No caso de barragem 
ou instalação em operação em que seja identificada comunidade na ZAS, deverá ser feita a 
descaracterização da estrutura, ou reassentamento da população e o resgate do patrimônio 
cultural, ou obras de reforço que garantam a estabilidade efetiva da estrutura, em decisão 
do poder público, ouvido o empreendedor e consideradas a anterioridade da barragem em 
relação à ocupação e a viabilidade técnico-financeira das alternativas. 

Na análise da segurança de barragens deve-se considerar, em atendimento 
a resolução nº 143/2012 do CNRH- conselho nacional de recursos hídricos: existência 
de documentação de projeto, incluindo o projeto “as-built”, estrutura organizacional e 
qualificação dos profissionais da equipe técnica de segurança da barragem; existência 
de planos de ações de emergência; procedimentos de inspeções de segurança e 
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monitoramento; relatórios de inspeção de segurança com análises e interpretação e 
regra operacional dos dispositivos de descarga da barragem; existência de população à 
jusante e potencial de perda de vidas humanas; existência de instalações de infraestrutura 
ou serviços e de equipamentos de serviços públicos essenciais. O tipo de barramento 
também é considerado importante, pois cada tipo de barragem possui sua susceptibilidade 
particular a determinadas falhas; em barragens de terra o maior problema concentra-se na 
possibilidade de ocorrência de erosão interna por piping, descrito por Ladeira (2007) como 
uma erosão interna regressiva (de jusante para montante). Nas barragens de concreto 
o mecanismo de falha pode se dar por meio de fissuras no concreto, na grande maioria 
das vezes resultado da má cura do concreto. Em obras que utilizam grandes volumes de 
concreto podem ocorrer fissurações ou trincas na estrutura devido a tensões no concreto 
geradas pelo gradiente de temperatura durante o processo de cura. 

O ISB- Índice de Segurança de Barragens considera 29 critérios para seu cálculo, 
segundo Zuffo (2005) os parâmetros considerados no seu cálculo são: altura da barragem, 
comprimento, tipo de barragem (material), tipo de vertedor, período de retorno da vazão de 
projeto dos vertedores, idade da barragem, presença de percolação/vazamento, presença 
de deformações e recalques, deterioração dos taludes, conservação das estruturas 
vertedoras, evidências de erosão a jusante. 

A proposta desse trabalho é propor o ISASB, Índice de Sustentabilidade Ambiental 
de Segurança de Barragens, um índice que seja analisado ao lado do ISB. E com a 
análise do ISB e do ISASB tem-se um critério para analisar a estabilidade do barramento, 
considerando também a questão ambiental e a bacia onde o barramento está incluído. 
Considerando os desastres ocorridos no Brasil e no mundo vários episódios e problemas 
são devidos a aspectos ambientais, como fortes chuvas, chuvas longas à montante do 
barramento. Assim um índice que considere tais parâmetros é importante. 


METODOLOGIA 


O ISB- Índice de Segurança de Barragens é dado pela fórmula dada na equação 1: 
n 
i=1 (1) 
Onde: qi: nota do i-ésimo critério, um número entre 1 e 100, obtido na respectiva 
curva da função de valor, em função do conceito obtido no critério; wi: peso correspondente 
ao i-ésimo critério, um número entre O e 1, atribuído em função da sua importância para 
a conformação global de segurança. Os pesos (wi) padronizados foram estabelecidos a 
partir de opiniões técnicas de diversos profissionais da área, obtidas a partir da resposta do 
questionário Aguiar et al. (2015), 
O índice de Sustentabilidade Ambiental de Segurança de Barragem- ISASB, 
proposto neste trabalho, considera critérios ambientais para seu cálculo e a utilização deste 
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índice leva em consideração fatores ambientais da bacia onde a barragem está incluída. 


A segurança de barragens é importante considerando os imensos desastres que já 


ocorreram no mundo e no Brasil. Entre estes desastres temos os abaixo relacionados, a 


maioria dos pesquisados com perdas materiais e humanas, tabela 1 e figura 2: 


ANO LOCAL NOME MORTOS/DESABRIGADOS/ CAUSA 
CONSEQUÊNCIAS 
1864 Inglaterra Barragem 250 mortos Fissura-percolação 
DaleDyke 
1889 USA- Barragem 2.209 mortos. Após chuvas intensasFoi o 
Pensilvânia-  SouthFork ondas de 12 metros maior desastre da história. 
Johnstown Rio Little chegaram aJohnstown, 20 milhões toneladas de água 
Conemaugh situada a 50 km da atingiram a cidade. (= toda a 
barragem, em 1Ominutos. água do rio mississipi) 
1923 tália- Barragem de 4,5 milhões de lama Físicas-Pressão 
Bérgamo Gleno atingiram Alterações feitas na Barragem 
as cidades 
1928 USA- Barragem 450 mortos Instabilidade Geológica do 
Califónia SaintFrancis 47 milhões litros de água Cânion. Estrutura “gravidade 
inundouo Vale de Santa arqueada”. Modificou-se a 
Clarita altura (mais 3 metros) para 
aumentar a capacidade., sem 
compensação na largura. 
1959 França- Barragem 423 mortos Separação da rocha da sua 
Frejus Malpasset 50 milhões de metros fundação, pela elevação da 
cúbicos de água a70 km/ água na rocha de fundação. 
hora. 
1960 BRASIL- Barragem de 1.000 mortos. Durante a construção. 620 de 
Ceará Orós comprimento e 54metros de 
altura. 
Fortes chuvas causaram 
o rompimento antes da 
conclusão. 
1975 China Bangiro e 26 mil mortos Período chuvoso de 24 hs e 
Shimantan 6 milhões casas destruídas - máconstrução 
145 mortos fome e 
epidemias 
1,7 bilhões de metros 
cúbicos de águaperdidos. 
1976 USA-Idah- Barragem de Barragem de Terra, com Ruptura por erosão laminar 
Teton Teton 123 me 356km?, gerou regressiva (piping). 
onda com 66.000 m?/s. 11 
vítimas e várias cidades. 
1979. Índia-Morbi Machchu-2 25 mil mortos Desmanchou-se após 


atingiu a cidade de Morbi, a 
5 km darepresa. 
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tempestade de 20 min com 9 
metros de altura 
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1980 


1985 


1992 


2000 


2008 


2010 


2015 


2019 


Espanha- 
Andraluzia 


Itália-Stava 


Srilanka 


USA 


Índia-Nepal- 
fronteira 


Hungria- 
Kolontar 


BRASIL-MG- 
Mariana 


BRASIL-MG- 
Brumadinho 


Aznalcollar 


Barragem Rio 
Stava 


Barragem de 
Kantale 


Martin County 
Barragem da 
Mineradora 
de Carvão 
Vegetal do 
sul dos USA 


Barragem 
Koshi-Nepal 


Barragem 
de Resíduos 
da Bauxita- 
Mineradora 
Mal, lama 
vermelha. 


Barragem 

de Fundão- 
Mineradora 
Samarco, Vale 
eBHP Biliton 


Barragem de 
rejeitos de 
mineração da 
Vale.Barragem 
B1- córrego do 
feijão. 


Tabela 1- Desastres com Barragens 


Avalanche de zinco, cobre, 
cádmio, sulfetos chegou aos 
rios Agro e Guadiamar. 

40 km de material tóxico 
atingiu o Parque Nacional 
de Dofiana, Patrimônioda 
humanidade da UNESCO. 
Contaminação do ambiente. 


268 mortos 

180 mil metros cúbicos de 
areia, água elama varreram 
a cidade de Stava. 


180 mortos 
8 mil afetados 


2 bilhões de lama de 
arsênico e mercúrioatingem 
os rios Tug Fork, Big Sandy 
e Ohio. 

Dissimou toda a vida 
aquática em poucos 
minutos. 


3 milhões desabrigados 
Centenas de mortos nos 
dois países. Áreas de 
agricultura soterradas 


milhão de metros cúbicos 
atingem ascidades de 
Kolontar e Devecser. 

4 mortos.Vários danos 
ambientais. Ondas de 3,5 
metros, barragem de 300 
metros de altura. 


62 milhões metros cúbicos 
de lama tóxica, atingem 
Bento Rodrigues emMinas 
Gerais. 

19 mortes. 

Várias perdas materiais. 
Perdas com Economia. 


12 milhões de metros 
cúbicos de rejeitosde 
minério. 

270 mortos e 11 
desaparecidos após 1ano 
do desastre. 
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Entupimento dos dutos de 
drenagem, com sedimentos, a 
água infiltrou-se no solo da 
represa. 


Ruptura da parede de 
contenção pelo excesso de 
tráfego pesado. 


Falta manutenção 


Processo Construtivo da 
Barragem. 


A Barragem rompeu. 
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Figura 2-Barragem de Saint Francis-USA. Antes e após a ruptura, em 192;. Barragem de Malpasset- 
França, antes e após a ruptura, 1959; e Distrito de Bento Rodrigues- Mariana, em Minas Gerais- 
Rompimento da Barragem de Fundão, em 2015. Fotos: Google(2021). 


As barragens são classificadas pelos agentes fiscalizadores, por categoria de risco, 
por dano potencial associado e pelo seu volume, com base em critérios gerais estabelecidos 
pelo Conselho Nacional de Recursos Hídricos (CNRH) (BRASIL, 2020). 

A classificação por categoria de risco em alto risco, médio e baixo risco será feita em 
função das caraterísticas técnicas, dos métodos construtivos, do estado de conservação e 
da idade do empreendimento e do atendimento ao plano de segurança de barragens, bem 
como de outros critérios definidos pelo órgão fiscalizador (BRASIL, 2020). 

A classificação por dano potencial associado à barramento em alto, médio e baixo 
risco será feita em função do potencial de perdas de vidas humanas e dos impactos 
econômicos, sociais e ambientais decorrentes da ruptura da barragem. O órgão fiscalizador 
deverá exigir do empreendedor a adoção de medidas que levem à redução da categoria de 
risco do barramento. 

São instrumentos da Política Nacional de Segurança de Barragens (BRASIL, 2010): 
O sistema de Classificação de barragens por categoria de risco e por dano potencial 
associado; O Plano de Segurança de Barragem, incluindo o PAE; o Sistema nacional de 
Informações sobre segurança de barragens (SNISB); o sistema nacional de informações 
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em meio ambiente (SINIMA); o cadastro técnico federal de atividades e instrumentos de 
defesa ambiental; o cadastro técnico federal de atividades potencialmente poluidoras ou 
utilizadoras de recursos ambientais; o relatório de segurança de barragens; o sistema 
nacional de informações sobre recursos hídricos (SNIRH); o monitoramento das barragens 
e dos recursos hídricos em sua área de influência e os guias de boas práticas em segurança 
de barragens, estes três últimos foram acrescentados pela Lei 14.006 de 2020 (BRASIL, 
2020). 

Aelaboração do PAE é obrigatória para todas as barragens classificadas como médio 
e alto dano potencial associado alto risco, a critério do órgão fiscalizador (BRASIL,2020). 
Mas independentemente da classificação quanto ao dano potencial ou ao risco, a elaboração 
do PAE é obrigatório para todas as barragens destinadas à acumulação ou a disposição de 
rejeitos de mineração (incluído na Lei 14.006 de 2020) BRASIL (2020). 


RESULTADOS 


A Matriz de categoria de risco e dano potencial associado para barragens, está na 
tabela 2: 


CATEGORIA DERISCO DANO POTENCIAL ASSOCIADO 
ALTO MÉDIO BAIXO 
ALTO A B [E 
MÉDIO A Cc D 
BAIXO A [6 E 


Tabela 2- Matriz risco e dano potencial associado (ANA, 2012) 


O relatório de segurança de barragens da ANA (2017) mostra que são mais de 7 
acidentes com barragens por ano no Brasil. Em 7 anos foram 52 acidentes registrados pela 
ANA- Agência Nacional de Águas. A tabela 3 mostra o registro anual desde 2011: 


2011 2012 2013 2014 2015 2016 2017 
ACIDENTE 4 1 - 5 4 6 4 
INCIDENTE 4 6 4 6 5 LE 10 


Tabela 3- Acidentes e Incidentes com Barragens de 2011 a 2017 (ANA,2017) 


Brunner (2014) apresenta uma lista de eventos ligados diretamente ao rompimento 
de Barragens: cheias, falhas de fundação, piping, falha na operação das comportas, 
problemas estruturais, escorregamento, abalos sísmicos, ruptura de uma barragem a 
montante, esvaziamento rápido do reservatório, sabotagem e remoção planejada do 
barramento. 

Para evitar que eventos como estes ocorram o PAE- Plano de Ação de Emergência 
estabelecerá as ações a serem executadas pelo empreendedor da barragem em caso 
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de situação de emergência, bem como identificará os agentes a serem notificados dessa 
ocorrência. 

Aguiar et al.(2015) adicionou os seguintes parâmetros ao cálculo do ISB: o volume 
do reservatório, a vazão de projetos dos vertedores; a conservação das estruturas de 
captação, a instrumentação e monitoramento dos registros; a existência de barragens em 
série no mesmo curso, de água (efeito cascata), tipo de turbina hidráulica, tipo de comporta 
do vertedor, maquinário de,operação das comportas, nível de automação e capacidade 
instalada ou em operação. 

Aguiar et al. (2015) sugeri a inclusão de outros parâmetros ao ISB : tipo de 
fundação, tipo de uso do barramento, alterações na bacia hidrográfica à frente às ações 
antrópicas, taxa de assoreamento do lago, existência de eclusas e onda de cheia a jusante. 
Todos estes não consideram a variável ambiental, que influencia muito a segurança das 
barragens, como parâmetros que consideram as áreas verdes ao redor dos reservatórios 
e a intensidade das chuvas à montante do barramento, que causam alterações na vazão 
do rio e carreamento de sedimentos que influenciam a qualidade das águas que chegam 
aos reservatórios, assim como o assoreamento existente no reservatório e que diminue a 
área útil deste. 

Assim este trabalho propõe a consideração dos seguintes parâmetros no cálculo 
deste novo índice, Índice de Sustentabilidade Ambiental de Segurança da Barragem- 
ISASB: 

Kas=O assoreamento existente na área do barramento: volume assoreado/volume 
morto da barragem, km3/kmS; 

Kav=A quantidade de área verde ao redor do barramento: km/km?; 

Kum=0O uso múltiplo da barragem: quantidade de usuários/km?; 

Kabd=Área degradada da Bacia de Drenagem: km2/km?; 

Kacvb= área de cobertura vegetal na bacia de drenagem: km?/km?; 

Kic= intensidade das chuvas na área da bacia, à montante da barragem: mm/hora. 


n 

ISASB = [ar (2) 
Onde: qi: nota do i-ésimo critério, um número entre 1 e 100, obtido na respectiva 
curva da função de valor, em função do conceito obtido no critério; wi: peso correspondente 
ao i-ésimo critério, um número entre O e 1, atribuído em função da sua importância para 
a conformação global de segurança. Os pesos (wi) padronizados serão estabelecidos a 
partir de opiniões técnicas de diversos profissionais da área. Os qi são critérios de 1 a 
100, relacionados aos parâmetros: Kas, Kav, Kum, Kabd, Kacvb e Kic. Wi são os pesos 

atribuídos a estes critérios. 
O Kas=0O assoreamento existente na área do barramento: volume assoreado/volume 
morto da barragem, km3/km3 é um parâmetro que permitirá avaliar a situação da barragem 
quanto ao assoreamento existente e como compromete a estrutura física da barragem; o 
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parâmetro Kav=a quantidade de área verde ao redor do barramento: km/km2 permiti avaliar 
como a barragem está sujeita a desmoronamentos de áreas próximas a ela, sendo que as 
áreas verdes permitem proteção do solo da área do entorno e evitam assoreamento do 
reservatório; o Kum=O uso múltiplo da barragem: quantidade de usuários/km? permite que 
seja analisada como é o aproveitamento e quantos dependem da barragem e contribuem 
para diminuir o seu nível, diminuindo os riscos de inundações; o Kabd=Área degradada 
da Bacia de Drenagem: km?/km? permite avaliar como está a área da bacia de drenagem 
e a situação de erosão e transporte para o reservatório provocando assoreamento; o 
parâmetro Kacvb= área de cobertura vegetal na bacia de drenagem: km2/km? permite 
avaliar como está a proteção de toda a bacia, onde a barragem está incluída e não apenas 
ao redor desta, como prevê a legislação; o Kic= intensidade das chuvas em 24 horas, na 
área da bacia, à montante da barragem: mmvhora, permite avaliar se a quantidade de água 
provoca inundações no reservatório, devido às variações do El Nifio e da La Nifia e outros 
fenômenos climáticos. 

Assim a proposta é que o Índice de Segurança de Barragens Global seja a soma: ISB 
+ ISASB. Considerando tanto os parâmetros físicos, estruturais da barragem (ISB) quanto 
os ambientais da bacia (ISASB), onde o barramento está inserido. Na figura 3 tem-se um 
exemplo de Barragem de abastecimento de água com cobertura vegetal protegida e com 
margens preservadas. A Água infiltra no solo devido a preservação da cobertura vegetal 
e florestal na Bacia e não provoca erosão e carregamento de partículas para a barragem, 
também tem-se a preservação da mata ciliar ao longo do rio e da barragem, ambos evitam o 
assoreamento desta. Existe um controle do uso e ocupação do solo próximo ao barramento 
o que reduz a interferência antrópica. A primeira barragem teve início em 1894, e foi 
concluída em 1918. Entre 1945 e 1950 outra barragem foi construída, com 30 a 50 metros 
de profundidade e fornece 20 milhões litros de água, sendo administrada pela CESAN, 
considerada em situação de estabilidade pela AGERH, Agência de Águas do ES (FRANÇA 
e COSTALONGA, 2013). Em 1965 virou reserva florestal e em 1991 Reserva Biológica, 
sendo as áreas ao redor desapropriadas em 1912, para a preservação da floresta. Assim o 
ISASB resulta em um bom índice. 


Figura 3- Barragem Duas Bocas em Cariacica- Reserva Biológica Duas Bocas. (IEMA, 2021) 
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CONCLUSÕES 


O Índice de Sustentabilidade Ambiental permite que os parâmetros ambientais sejam 
levados em consideração quando da análise dos riscos a desastres com barramentos, 
principalmente com núcleos urbanos situados à jusante destes. 

Ações de reflorestamento e recuperação da cobertura vegetal ao redor destes 
barramentos possibilita ampliar a capacidade de armazenamento e resulta em índices de 
segurança melhores, pois contribuem para a infiltração de água no solo, a proteção do solo 
contra a erosão e o assoreamento do barramento. 

Um desafio é criar critérios que considerem estes parâmetros no cálculo do ISASB 
(gi e wi), o que necessita da experiência dos especialistas da área e de pesquisas. Outro 
grande desafio é a Educação Ambiental das comunidades ao redor do Barramento e a 
Gestão Ambiental da bacia Hidrográfica onde ele está inserido. Uma sugestão é o trabalho 
com os Comitês de Bacia Hidrográfica e a ampliação da participação de usuários, poder 
público e sociedade civil. Outra sugestão é investir num moderno Sistema de Informações 
Georreferenciadas e disponíveis e numa rede de monitoramento eficiente. 
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CAPÍTULO 2 
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George Rodrigues de Sousa Araújo 
Serviço Geológico do Brasil, Recife, PE, 
Brasil 


André Luiz Andrade Simões 


Universidade Federal da Bahia, Salvador, 
BA, Brasil 


Rodrigo de Melo Porto 


Escola de Engenharia de São Carlos — 
Universidade de São Paulo, São Carlos, 
SP, Brasil 


RESUMO: A Rede Hidrometeorológica 
Nacional - RHN monitora as vazões dos 
principais rios brasileiros para subsidiar o 
planejamento da utilização dos recursos 
hídricos em todo território nacional e utiliza 
muito frequentemente o equipamento 
“Perfilador 


Corrente por efeito Doppler” — ADCP em 


denominado Acústico de 
suas operações hidrométricas. O principal 
objetivo deste trabalho foi comparar as 
vazões calculadas empregando o ADCP 
em modo estático e o princípio da máxima 
entropia, com as vazões medidas pelo 
ADCP em modo dinâmico. foram realizadas 
medições de vazão em vinte estações 
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situadas no estado da Bahia em Sergipe. 
O intervalo de confiança do parâmetro de 
entropia M foi estimado valendo-se dos 
registros históricos de velocidades pontuais 
medidas com molinete hidrométrico e 
os demais parâmetros da metodologia 
alternativa proposta 
ao perfil de máxima velocidade pontual 
por regressão não linear. Os resultados 


foram ajustados 


evidenciaram que em 79,4% das medições, 
os valores calculados de vazão divergiram 
em até +10% em relação aos medidos com 
o ADCP em modo dinâmico; o tempo médio 
de medição, em campo, decresceu 46,34%; 
ausência de limite físico na magnitude da 
vazão para aplicação da metodologia da 
máxima entropia no cálculo da descarga 
líquida em canais naturais. 
PALAVRAS-CHAVE: Medição de vazão; 
Entropia; ADCP. 


ABSTRACT: The National 
Network - RHN 
monitors the discharge of the main brazilian 
rivers to subsidize the planning of the use 
of water resources throughout the national 


Hydrometeorological 


territory and very often uses the equipment 
called “Acoustic Doppler Current Profile” - 
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ADCP in its hydrometric operations. The main objective of this work was to compare the 
calculated discharge rates using ADCP in static mode and the principle of maximum entropy 
with the discharge measured by ADCP in dynamic mode. 34 discharges were measured in the 
interval 6.138 mº/s - 1,326.8 m/s, distributed in 20 RHN stations in the state of Bahia and 1 in 
Sergipe. The confidence interval ofthe entropy parameter M was estimated using the historical 
records of point velocities measured with an analog current meter, the other parameters of the 
proposed alternative methodology were adjusted to the profile of maximum point velocity by 
non-linear regression. The results showed that: in 79.4% of the measurements, the calculated 
values of flow diverged up to + 10% in relation to those measured with ADCP in dynamic 
mode; the mean measurement time, in the field, decreased by 46.34%; absence of physical 
limit on flow magnitude for application of the maximum entropy methodology in the calculation 
of discharge in rivers. 

KEYWORDS: Discharge Measurement, Entropy, ADCP. 


INTRODUCTION 


Escoamentos turbulentos de fluidos newtonianos são tridimensionais, dependentes 
do tempo e formados por estruturas internas que compõem uma variedade de escalas, 
desde as grandes escalas até as microescalas de Kolmogorov. A sua representação física 
— matemática está consolidada desde 1845, quando Stokes publicou o seu artigo com 
a versão final das equações hoje conhecidas como equações de Navier-Stokes. Para a 
turbulência, essas equações em conjunto com a equação de conservação de massa não 
possuem solução analítica conhecida, condição que estimulou a proposição de modelos 
para a turbulência, cujos níveis vão desde a modelagem direta das covariâncias que 
compõem o tensor de Reynolds às representações mais restritas apenas às distribuições 
de velocidades. 

Os modelos de distribuição de velocidades em escoamentos em condutos forçados 
e condutos livres têm origem em análise estatística de dados experimentais, considerações 
fenomenológicas e análise dimensional, a exemplo das equações logarítmicas e leis de 
potência. A formulação logarítmica, amplamente utilizada nos estudos de Prandtl, von 
Karman e Nikuradse representa uma importante contribuição ao tema, mas não descreve 
a distribuição próxima a parede. Nesse contexto, Chiu (1987) propôs uma distribuição sem 
esta restrição, tendo como base a teoria da informação de Shannon (1948). O emprego da 
formulação proposta por Chiu (1987, 1988), assim como ocorre com os outros modelos, 
requer a determinação de um parâmetro de ajuste. O conhecimento desse parâmetro para 
uma determinada estação fluviométrica proporciona mais uma alternativa à determinação 
da vazão, além de conduzir a uma metodologia cujo tempo de medição é relativamente 
reduzido. O presente estudo teve como objetivo a determinação do referido parâmetro 
da formulação baseada na entropia da informação, para dados experimentais obtidos em 
21 estações fluviométricas pertencentes à Rede Hidrometeorológica Nacional - RHN: 20 
localizadas no estado da Bahia e 1 em Sergipe. Também foi objetivo deste trabalho o 
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cálculo da vazão com a formulação baseada na entropia e sua comparação com vazões 
medidas utilizando ADCP e a análise do tempo de execução de cada método, tendo em 
vista a proposição de um método alternativo de medição de vazão com tempo reduzido em 
relação ao convencional, baseado no uso do ADCP. 


MATERIAL AND METHOD 


A formulação empregada neste trabalho para determinação modelagem da 
distribuição de velocidades e cálculo de grandezas relacionadas foi apresentada por Chiu 
(1988). Sua dedução e aspectos conceituais podem também são discutidos em Moraes 
(2010), Araújo (2017) e Martins e Poleto (2017). 

Em conformidade com Chiu (1991), sendo o eixo y a vertical onde se encontra o 
ponto de máxima velocidade de escoamento na seção, u 
é expressa por: 


a distribuição de velocidades 


max? 


NE ma |n (1 + (eM — Det) (1) 


em que, D é a altura de escoamento na seção transversal, h é a distância vertical 
desde a superfície livre até o ponto de máxima velocidade e M é o parâmetro de ajuste, cuja 
definição está trelada à velocidade máxima e a um multiplicador de Lagrange empregado 
na dedução. 

Uma segunda formulação resultante do trabalho de Chiu (1988) empregada neste 
trabalho é a equação 2, que relaciona a velocidade média, a máxima e M. Esta equação é 
especialmente útil para determinação da vazão com o uso da velocidade média. 


M 1 


= -=6 (2) 


Umáx = (eM-1) M 


El 


Seleção das estações fluviométricas 


Foram selecionadas para esta pesquisa 21 estações fluviométricas pertencentes à 
RHN: 20 localizadas no estado da Bahia e 1 em Sergipe. Estes postos fluviométricos foram 
selecionados por possuírem profundidade média compatível com a utilização do Perfilador 
Acústico de Corrente por efeito Doppler — ADCP nos períodos de vazões máximas e 
mínimas. As estações selecionadas estão elencadas na Figura 1 e na Tabela 1. 
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adaptado de Araújo (2017). 
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Figura 1. Localização das estações fluviométricas selecionadas para medição de vazão com ADCP, 


Estação Nome Rio a E 
45770000 Arrojado Arrojado 5540 
46550000 Barreiras Grande 24400 
45480000 A São Francisco 271000 
46902000 Boqueirão Grande 65900 
45590000 Correntina Correntina 3900 
50191000 e Vaza Barris 15740 
46870000 fazenda Preto 22000 
Porto Limpo 

46035000 Gameleira São Francisco 309000 
46830000 Ibipetuba Preto 17900 
46150000 Ibotirama São Francisco 323000 
48020000 Juazeiro São Francisco 510800 
45740001 Mocambo Do Meio 7950 
46360000 Morpará São Francisco 345000 
46105000 Paratinga São Francisco 314000 
45960001 Porto Novo Corrente 29700 
45910001 Santa Maria Corrente 28200 
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46650000 Taguá Grande 36000 


50795000 Tiririca Pojuca 4860 
51795000 Valença Una 1110 
46780000 Vereda Preto 1900 
45580000 Veredão Corrente 1472 


Tabela 1. Estações selecionadas para pesquisa. 


Intervalo de confiança de Ze M 


Nas estações da RHN selecionadas para esta pesquisa, foram analisadas as 
velocidades pontuais, medidas com molinete hidrométrico, em formulários originais 
de medição de vazão a partir do ano 2001, digitalizados e disponíveis no Sistema de 
Visualização de Ficha de Campo da ANA, disponíveis em ANA (2016). 

As etapas de cálculo para estimativa de M foram: 


* Registro das velocidades média e máxima em todos os formulários de uma 
mesma estação, computando-se o valor de 9 em cada medição de acordo com 
a Equação 2; 


*-— Admitindo-se que os valores de 9 tendem a seguir uma distribuição normal de 
probabilidades, determinou-se o intervalo de variação de com base no Intervalo 
de Confiança-lC (95%) sobre a média (9) em cada estação, de acordo com a 
expressão: 


QD — 4 S0' Fo) ml 3 
M(6 — 1,96 o? <M< M(6+ 1,96 a , (3) 


em que, para cada estação: 9=média amostral dos valores de 9, o =Desvio padrão 
amostral de 9 e n,=tamanho da amostra de o. 


«Foram calculados os valores de M(9) com o método de Newton-Raphson. 


Medição de vazão 


Para medir a vazão e caracterizar o perfil vertical de máxima velocidade das estações 
fluviométricas selecionadas, foi empregado um equipamento ADCP com quatro feixes e de 
dupla frequência M9 (1 e 3 MHz). 

A vazão média e a área da seção foram obtidas empregando o ADCP em modo 
dinâmico e seguindo as recomendações de ANA (2014): 


e Travessias aos pares, com igual número de travessias em ambos os sentidos 
da seção; 


* Tempo de medição de no mínimo 720 segundos; 
* Velocidade do barco não superior ao dobro da água; 


* | Coleta de 10 verticais nas margens iniciais; 
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*. Trajeto de medição perpendicular à seção e alinhado; 


* Verificação da ocorrência do efeito de fundo móvel. 


Considerando a necessidade deidentificar avertical de máxima velocidade e questões 
de ordem prática e econômica, dividiu-se a seção transversal em subseções, semelhante à 
medição com molinete, para execução do perfilamento vertical de velocidades com o ADCP 
em modo estacionário, de acordo com as sugestões de Collischonn & Dornelles (2013) 
contidas na Tabela 2. 


Largura do rio Número de 
(m) Verticais 

<=30.00 8 
30.00-50.00 10 
50.00-80.00 13 
80.00-150.00 14 
150.00-250.00 18 
>= 250.00 20 


Tabela 2. Distância de perfilamento entre as verticais, adaptado de Collischonn & Dornelles (2013). 


O tempo de perfilamento por vertical foi de 60 segundos. Tanto o tempo necessário 
para executar as travessias da medição em modo dinâmico quanto o tempo total de 
perfilamento em modo estacionário foram registrados pelo ADCP. 


Ajuste numérico 


Selecionou-se a vertical cuja velocidade média temporal é máxima durante a 
medição com ADCP em modo estacionário. A Equação 4 expressa, matematicamente, o 
ajuste da Equação 1 ao conjunto de pontos (u, y) da vertical de máxima velocidade média 
temporal na seção: 


minimizar (Nº. ,(u; — u)?), (4) 


em que, 
Um ” yo (1-5 ' 
u = u(yi) = “kn (1 + (e = o D-h ) » (5) 


e sujeito à restrição da Equação 3. 

Não houve restrições de variação em relação âu, , eh. 

Foram calculados também os coeficientes de determinação (R?) entre as velocidades 
medidas (ADCP em modo estacionário) e as calculadas com a Equação 1. 

Após o ajuste de mínimos quadrados obteve-se o valor de M que melhor descreve 
o comportamento da distribuição de velocidades na seção; a Equação 2 foi então utilizada 
para calcular g(M), possibilitando o cálculo da vazão com 
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Qcar = A*U=Uumtb, (6) 
em que: Q.. =Vazão na seção calculada pelo método da máxima entropia e A= área 


da seção transversal. 


Calc 


Análises estatísticas 


Calcularam-se erros relativos percentuais entre as vazões medidas com ADCP e as 
calculadas com a Equação 6. 


. si: (Q alc—Q e ) 
Difg = Sete Sted: «100, (7) 


em que: Dif, = erro relativo entre as vazões, Q 
da entropia, Ques 
A divergência entre os tempos de medição está de acordo com a Equação 8. 


= vazão calculada pela metodologia 


Calc 


= vazão medida com o ADCP em modo dinâmico. 


Difr = Sino) 100), (8) 

em que: Dif = divergência relativa entre os tempos de medição em campo, T,,=tempo 
de medição empregando a metodologia da entropia com ADCP em modo estacionário, 
Tom=iempo de medição da vazão com ADCP em modo dinâmico. 

Com propósito de verificar uma potencial diferença significativa entre: a média (pg) 
e a variância (02) das vazões de amostras distintas (ADCP estacionário/entropia e ADCP 
dinâmico), e o tempo de medição obtido por diferentes metodologias, aplicaram-se análises 
estatísticas: 


* | Regressão linear entre as vazões medidas (Q,,.,) e calculadas (Q,...); 


* Correlação linear entre a diferença relativa percentual de vazões (Dif,) e as 


variáveis Quea € Dify; 


* Teste de significância t-student, com dados pareados, para igualdade da média 


das vazões das variáveis Q.. € Quea: 


Calc 
* Teste de significância F-Snedecor, com dados pareados, para igualdade da va- 


riância das vazões das variáveis Q,.. e Q 


Cale Med” 


* Teste de significância t-student (pareado) para igualdade do tempo médio de 
execução das medições, Tp, € Tem 


RESULTADOS E DISCUSSÕES 


As Tabelas 3 e 4 resumem os resultados da análise da série de velocidades máximas 
e médias medidas com molinete nas estações com propósito de determinar o intervalo para 
o parâmetro M no ajuste numérico. 
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Estação Reno 


histórica 


no uv(o) (2) 


Mm 
5 (4%) 


45770000 2001-2013 35  0,6449 0,0481 
46550000 2001-2013 34 0,7121 0,0252 
45480000 2001-2013 34 0,6643 0,0767 
46902000 2001-2013 31 0,7378 0,0443 
45590000 2001-2013 31 0,5740 0,0380 
50191000 2001-2013 50 0,6755 0,0914 
46870000 2001-2013 40 0,6297 0,0560 
46035000 2001-2013 35 0,7165 0,0537 
46830000 2001-2016 47 0,6844 0,0439 
46150000 2001-2013 36 0,7458 0,0482 
48020000 2002-2013 25 0,7209 0,0363 
45740001 2001-2013 34 0,5899 0,0612 
46360000 2001-2013 28 0,6719 0,0643 
46105000 2001-2013 33 0,7340 0,0471 
45960001 2001-2013 37 0,7918 0,0338 
45910001 2001-2013 30 0,7068 0,0212 
46650000 2001-2013 33 0,7924 0,0216 
50795000 2001-2012 33 0,5919 0,1044 
51795000 2001-2015 48 0,6092 0,0996 
46780000 2002-2016 38 0,6497 0,0264 
45580000 2001-2013 37 0,6331 0,0331 


7,46 
3,53 
7,67 
6,00 
6,63 
13,54 
8,89 
5,37 
6,41 
4,82 
5,03 
10,38 
9,57 
6,41 
4,26 
3,00 
2,73 
17,63 
16,35 
4,06 
5,28 


Tabela 3. Estatística descritiva dos valores de 9 por estação, elaborado com dados da ANA (2016). 


Estação D/1C-95% IC-M(9) 
45770000 | 0,628-0,661 1,608-2,067 
46550000  0,703-0,720 2,730-3,030 
45480000  0,638-0,690 1,737-2,520 
46902000  0,721-0,753 1,737-2,520 
45590000  0,560-0,587 1,737-2,520 
50191000 | 0,650-0,700 3,052-3,672 
46870000 | 0,611-0,647 3,052-3,672 
46035000  0,698-0,734 0,730-1,071 
46830000  0,671-0,697 0,730-1,075 
46150000  0,729-0,761 -0,684-5,234 
48020000  0,706-0,735 1,386-1,870 
45740001 0,568-0,610 2,644-3,287 
46360000  0,647-0,696 2,640-3,287 
46105000  0,701-0,750 2,644-3,287 
45960001 0,780-0,802 2,222-2,627 
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45910001  0,699-0,714  3,195-3,847 
46650000  0,784-0,799  3,195-3,847 
50795000  0,556-0,627  3,195-3,847 
51795000  0,580-0,637 | 1,618-4,636 
46780000  0,641-0,658  0,836-1,972 
45580000  0,622-0,644 | 1,871-2,609 


Tabela 4. Intervalo de variação de M para o ajuste numérico, elaborado com dados da ANA (2016). 


Comparando o valor da razão (9) média entre a velocidade média e a máxima em 
cada estação nota-se que esse número varia entre 0,5740-0,7924 e o desvio padrão entre 
0,0212-0,0996. Em 33% das estações o Coeficiente de Variação (!!), CV, foi inferior a 5%, 
em 52% das amostras esse mesmo coeficiente oscilou entre 5%-11% e nas 15% estações 
remanescentes o CV oscilou no intervalo 11%-18%. Esses valores podem ser considerados 
de baixa dispersão, pois em 90,47% as amostras apresentam valores de CV inferior a 15%. 
Após o cálculo do intervalo de variação de M empregando-se a Equação 3, os resultados 
dos ajustes numéricos da Equação 1 aos perfis de máxima velocidade pontual em cada 
estação podem ser visualizados na Figura 2 e na Tabela 5. 
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Figura 2. h= distância ao leito, D= profundidade da vertical de máxima velocidade pontual, V= 
velocidade calculada pela Equação 1, V MAX = velocidade máxima pontual medida com ADCP em 
modo estacionário. 


Observando a Figura 2, nota-se há diferentes quantidades de pontos medidos por 
estação, isso ocorre em função do ADCP discretizar as verticais em número de pontos 
distintos em cada uma delas, em função de diferentes profundidades médias, velocidades 
e sedimentos em suspenção. Perfis naturais aproximadamente monotônicos só foram 
encontrados na estação de 46902000-09/06/2016 e em 4559000-(26/10/2016), evidenciando 
a pequena probabilidade da ocorrência deste perfil naturalmente. Em 13 perfis analíticos, 
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a Equação 1 calculou a posição da velocidade máxima na superfície livre (nas estações 
46550000-02/06/2016, 45480000-14/06/2016, 46902000-04/05/2017). Nestes cenários, 
nota-se: o valor da variável h melhor ajustado aos dados é negativo; o valor da velocidade 


máxima (melhor ajustada) é superior ao valor da velocidade máxima real na seção; a 


velocidade máxima real é aquela calculada aplicando a Equação 1 à superfície livre no 


eixo y. A Equação 1 também foi matematicamente convergente ao valor nulo de velocidade 


nas regiões próximas ao leito. Ainda que a estação 51795000 (Valença) esteja localizada 


em região costeira, não foram encontrados valores negativos de velocidades medidas nas 


estações, ou seja, velocidades escoando no sentido de jusante para montante da seção. 


Analisando conjuntamente todos os perfis adimensionalizados de velocidades máximas 


medidas e calculadas da Figura 2, percebe-se que a Equação 1 pode ser considerada uma 


ótima ferramenta para descrever os perfis de velocidade verticais nas seções estudadas. 


Melhor ajuste 


A B 
Estação Ds Qua Ele Dif (m) dm) D Ton Ten Dif, V MÁX Vai pm Voa » h 
ms ms a m min min (%) m/s m/s m/s m 
ASTTOO0O | 20/10/2016 32,54 33,78 +38 338 38 1,18 23 16 -80,4 1,557 1,536 0,983 1,563 1,903 0,339 
46550000 02/06/2016 38,81 40,00 +31 57,2 38 1,80 22 17 -22,7 0,984 0,969 0,940 0,996 3030 1,010 
45480000 14/06/2016 411,00 448,68 +92 8980 301 351 57 27 -52,6 0,725 0,723 0,913 0,759 2520 2,586 
45480000 05/12/2016 123930 147,34 TA 17682 524 431 39 27 -30,8 1,008 0,981 0,948 0,981 2,068 1,238 
45480000 15/05/2017 417,41 392,82 59 TT 314 356 26 26 0,0 0,7767 0,731 0,559 07314 2,520 0,798 
46902000 09/06/2016 136,00 126,68 69 201,0 132 2,04 32 19 -40,6 0,848 0,836 0,970 0,836 3,672 0,166 
46902000 04/05/2017 | 172,08 171,15 05 2451 136 243 15 15 0,0 0,973 0,965 0,985 0,967 3073 — -0,346 
45590000 02/05/2016 20,77 17,38 163 253 27 1,05 25 1 -56,0 1,195 1,187 0,994 1,187 0,946 0,364 
45590000 26/10/2016 21,01 19,49 7,3 25,6 31 0,92 35 13 -62,9 1,334 1,294 0,901 1,294 1,075 0,250 
50191000 05/06/2017 6,13 6,07 41 9,8 19 071 17 20 +17,6 0,907 0,877 0,817 0,877 2,687 0,249 
46870000 09/06/2016 62,40 53,10 149 1410 49 383 35 13 -62,9 0,621 0,582 0,803 0,582 1,870 1,008 
46035000 15/06/2016 507,58 499,47 17 14390 394 472 2 33 -70,5 0,481 0,473 0,936 0,477 3280 1,326 
46035000 06/12/2016 132680 138141 +41 20680 405 5,83 35 24 -314 0,955 0,957 0,881 0,994 2640 — -3,753 
46035000 10/05/2017 | 562,97 548,97 -25 10780 358 5,63 27 27 0,0 0,722 0,693 0,703 0,693 3,280 0,585 
46830000 29/10/2016 69,97 72,73 +39 830 50 1,99 13 19 +46,2 1,301 1,305 0,958 1,306 2,222 0,150 
46150000 17/06/2016 507,60 526,55 +37 12070 335 547 99 34 -65,7 0,556 0,577 0,977 0,800 3735  -32,851 
46150000 08/12/2016 124428 126960 +20 18702 527 7,20 36 27 -25,0 0,920 0,891 0,494 0,897 3840 1,938 
46150000 0905/2017 | 549,63 543,20 12º 12214 332 581 27 28 +3,7 0,599 0,609 0,904 0,736 3,201  -15,843 
48020000 01/06/2017 | 684,91 70688 +32 21270 685 303 50 27 -46,0 0,465 0,4520 0,890  0,4521 3,299 — -0,100 
45740001 | 25/10/2016 19,22 18,83 -20 34,6 49 09 30 14 -53,3 0,944 0,954 0,991 0,954 0,837 0,171 
46360000 20/06/2016 525,00 455,14 133 8780 366 350 87 33 -62,1 0,780 0,745 0,757 0,745 2,609 0,947 
46360000 09/12/2016 124212 112200 -97 15162 428 374 30 24 -20,0 1,104 1,063 0,906 1,063 2609 — -0,094 
46360000 12/05/2017 | 544,88 488,70 403 9114 393 306 30 22 -26,7 0,787 0,770 0,918 0,770 2,609 0,645 
46105000 16/06/2016 511,00 716,38  +40,2 9920 596 360 142 47 -66,9 1,041 1,030 0,898 1,030 2,690 0,295 
46105000 07/12/2016  1263,14 126811 +04 18797 610 312 45 38 -15,6 0,945 0,908 0,815 0,908 3,458 0,429 
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46105000 


45960001 


45910001 


45910001 


46650000 


50795000 


51795000 


46780000 


45580000 


11/05/2017 555,59 596,98 +74 1141,9 605 2,46 49 31 -36,7 0,765 0,7457 0,981 0,7673 2,690 


22/04/2016 105,83 111,82 +5,7 234,9 66 4,02 36 16 -55,6 0,635 0,610 0,719 0,610 4,287 
18/04/2016 110,86 111,46 +0,5 179,9 99 2,23 39 18 -53,8 0,892 0,886 0,970 0,886 2,650 
02/06/2017 115,60 129,41 +11,9 167,5 88 3,08 30 19 -36,7 0,936 1,105 0,968 5,091 2,650 
01/06/2016 85,85 90,57 +5,5 109,6 87 1,89 46 17 -63,0 1,042 1,053 0,949 1,053 4,390 
07/06/2017 115,60 129,41 +3,1 34,6 44 1,01 13 17 +30,8 0,318 0,313 0,958 0,3201 4,130 
05/09/2016 24,31 25,50 +4,9 110,5 36 3,58 44 14 -68,2 0,376 0,362 0,827 0,362 1,730 
27/10/2016 22,98 26,96 +17,3 40,7 29 2,23 14 12 -14,3 1,077 1,026 0,909 1,026 1,831 
29/04/2016 13,24 13,76 +3,9 16,2 13 1,50 30 n -63,3 1,251 1,361 0,927 1,416 1,522 

n 34 34 34 34 

V 401,027 402,320 40,8 22,2 

o 428,778 425,018 28,6 85 


-1,333 


0,603 


0,268 


-135,41 


-0,001 


-0,580 


1,214 


0,357 


-0,777 


B= largura superficial, Va = Velocidade máxima calculada pela Equação 17 na seção, V 
=velocidade máxima resultante do ajuste numérico, R2= Coeficiente de determinação entre velocidades 
medidas e as calculadas pela equação 1. 


Tabela 5. Resumo dos principais parâmetros do ajuste da Equação 1 aos perfis de máxima velocidade. 


De acordo com a variável V. MÁX da Tabela 5, percebe-se uma distribuição das 
velocidades máximas medidas no intervalo 0,318-1,557 m/s com prevalência de velocidades 
de magnitude 0,5-1,0 m/s representando mais de 55% da amostra; considerando o intervalo 
0,5-1,5 m/s, têm-se 83% da amostra. Estes números estão consoantes com os intervalos 
de dados de velocidades máximas obtidas nos experimentos de Minei (1999), Chen et al. 
(2013) e Corato et al. (2014). 

Os valores de M melhor ajustados aos dados reais estão no intervalo [0,837-4,39]. 
Comparando-se os valores de M calculados neste trabalho com os obtidos por diferentes 
pesquisadores percebe-se que aqueles valores são compatíveis com estes, como 
demonstrado por Minei (1999), Chen (2013), Choo et al. (2013), Farina et al. (2014), cujos 
valores de M calculados estão no intervalo [1,4-73,63], [0,139-0,206], [4,558-7,71], [1,26- 
4,08], respectivamente. 

Considerando as dezesseis medições realizadas no rio São Francisco como uma 
só amostra, nota-se que o valor do parâmetro da entropia (M) melhor ajustado aos perfis 
de velocidade permaneceu no intervalo [2,609-3,840]. Esses dados sugerem que para 
estações de mesma bacia hidrográfica os valores do parâmetro M podem variar em um 
especiro específico e característico da própria região hidrográfica. 

Os coeficientes de determinação (R?), resultantes do ajuste entre os dados de 
velocidade medidas (ADCP estacionado) e as calculadas pela Equação 1, variaram no 
intervalo [0,494-0,994]. Analisando a Figura 2 e a Tabela 5, percebe-se que em 64% (22) 
das estações analisadas o valor de Rº superou 0,9, resultando em um ótimo ajuste entre as 
velocidades medidas e calculadas; em 20,5% (7) das medições o valor de R? permaneceu 
no intervalo [0,7-0,8], consideradas de razoável ajuste; e nas estações restantes (5, ou 
14,7%) os ajustes foram considerados ruins. 
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Conforme a Figura 3, as vazões medidas e calculadas estão fortemente alinhadas e 
correlacionadas (R?=0,9864), principalmente para vazões inferiores a 200 m/s. A Figura 3 
sinaliza que as vazões calculadas pela metodologia alternativa da entropia podem computar 
vazões em substituição à metodologia do ADCP dinâmico. 


Qmed X Qcaic 


y=0,9931x 


3 
Rº = 0,9864 Quea (m?/s) 


Figura 3. Regressão linear entre vazões medidas e calculadas. 


Analisando a Figura 4, percebe-se que essas diferenças relativas entre as vazões se 
distribuem aleatoriamente no intervalo +20%, sendo a única exceção a estação 46105000 
em 16/06/2016, cuja diferença atingiu 40,2%. Analisando a distribuição das diferenças 
relativas das vazões em função das vazões medidas, percebe-se que não houve correlação 
(R2=0,004) entre essas variáveis. Portanto, para as vazões pesquisadas (6,138-1326,800 
mê/s), nota-se concordância com Diniz (2013), que também não detectou limitações para 
vazão em relação ao uso da metodologia da entropia para determinação da vazão. 


Quea X Difa (%) 
50 
40 
30 
20 + 
o 0% 
5 o É > 
a o 
-20 
-30 
-40 
-50 
o 200 400 600 800 1000 1200 1400 
y = -0,0016x + 1,5923 3 
m/s 
R? = 0,0046 Qua (m2/5) 


Figura 4. Correlação linear entre as variáveis Qu. € Difa- 
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De acordo com a Figura 5, a metodologia alternativa proposta reduziu o tempo 
de execução em 88,3% das medições, mas em 11,7% da amostra acarretou aumento 
neste tempo. Nota-se na Figura 5 que houve fraca correlação linear (R2=0,002) entre as 
variáveis Dif, e Dif,, sugerindo que a metodologia da entropia fornece dados de vazão com 
divergência aceitável (+10%) em relação à metodologia tradicional em 79,41% das estações 
amostradas independente da redução proporcional do tempo de medição, indicando um 
ganho de eficiência em quaisquer magnitudes de vazões a serem medidas. 


Dif, X Difo 
50 
40 
30 
20 
br 
O 10 
-20 2 
-30 
-40 
-50 
-80 -60 -40 -20 0 20 40 60 80 
y=-0,0157x + 0,4458 Dif, (%) 
R? = 0,0022 


Figura 5. Correlação linear entre as variáveis Dif, e Difa. 


Considerando uma amostra com 33 graus de liberdade (gl), a Tabela 6 indicou 
uma estatística do teste (t-calc=0,101) com probabilidade de significância (p-valor=0,918) 
muito superior ao nível de significância (a =0,05) adotado. Assim, a hipótese nula não deve 
ser rejeitada, e conclui-se que não há diferença significativa entre as médias das vazões 
medidas e calculadas. 


Hipótese nula (HO) 
Vazão média medida = Vazão média calculada 


Hipótese alternativa (H1) 
Vazão média medida * Vazão média calculada 


u (Difo) o (Dify) n gl tcalc a p-valor 
0,868 49,868 34 33 0,101 0,05 0,918 


Tabela 6. Teste de significância t-student (média das vazões) 


A estatística do teste F, F-calc=1,015, conduz a um p-valor (0,965) superior ao valor 
de a=0,05, indicando que a hipótese nula (igualdade de variâncias nas amostras) não pode 
ser rejeitada. Este teste sinaliza que, na amostra pesquisada, não há evidência suficiente 
para concluir que a metodologia alternativa da entropia estime vazões com variância 
diferente em relação ao ADCP dinâmico. 
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Os resultados apresentados na Tabela 8 mostram que ao valor da estatística do teste, 
t-calc=-4,704, está associado um p-valor muito próximo à nulidade (0,000022) e inferior ao 
nível de significância. Assim, conclui-se que a hipótese nula deve ser rejeitada em favor da 
hipótese alternativa. Há evidência estatística suficiente para concluir que o tempo médio de 
execução de uma medição de vazão utilizando o ADCP em modo estacionário (associado à 
metodologia da entropia) é menor em relação ao tempo dispendido com o ADCP em modo 


dinâmico. 
Hipótese nula (HO) 
Tempo-medição entropia = Tempo-medição dinâmica 
Hipótese alternativa (H1) 
Tempo-medição entropia < Tempo-medição dinâmica 
uv (Dif)  o(Dif) nm gl tcalc a p-valor 
-18,647 23,110 34 33 -4,704 0,05 0,000022 
Tabela 8. Teste de significância t-student (tempo de medição) 
CONCLUSÕES 


Os valores do parâmetro M da distribuição de velocidades foram ajustados com 
excelentes coeficientes de correlação e boa aderência entre teoria e experimentação. 
A aplicação da formulação demonstrou haver concordância com as vazões medidas 
empregando o ADCP. O uso da metodologia alternativa proposta para a determinação da 
vazão foi capaz de calcular a vazão em 79,4% dos casos, com divergência máxima de 
+10% em relação à metodologia convencional, e com redução média de 46,34% do tempo 
de medição. 

A metodologia da entropia com ADCP estacionário pode ser considerada uma 
alternativa para execução de medições regulares de vazão em redes hidrométricas, 
principalmente pela potencial redução do tempo de execução das medições em relação à 
metodologia com ADCP empregado em modo dinâmico. 
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RESUMO: Este artigo tem como principal 
objetivo apresentar as taxas de acesso 
à água tratada no mundo e no continente 
mostrar a 
águas 
meteóricas e incentivar seu aproveitamento 


americano, assim como 


importância da captação das 


a fim de sanar necessidades básicas de 
grande parte da população que não têm 
acesso à água de qualidade. Para isso, 
foi utilizado o método de estudo teórico 
documental com base no levantamento 
de dados disponibilzados por grandes 
organizações mundiais. Os resultados 
apresentaram concordância do uso de 
águas meteóricas, e em especial as águas 
pluviais, para fins potáveis e não potáveis 
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dependendo do uso, após tratamento 
físico-químico. Por fim, conclui-se que o 
investimento no aproveitamento de águas 
meteóricas pode auxiliar no acesso à água 
tratada para pessoas que não têm acesso 
à água de qualidade, bem como aliviar a 
demanda de água dos sistemas públicos de 
abastecimento de água trazendo economia 
de custo para a população beneficiada. 
PALAVRAS-CHAVE: Captação de água 
de chuva. Qualidade de água. Potencial de 
economia de água. 


METEORIC WATER AS A 
CONTRIBUTION TOWARD THE 
UNIVERSALIZATION OF ACCESS 
TO DRINKING WATER AND THE 
RESILIENCE TO THE EFFECTS OF 
CLIMATE CHANGE 
ABSTRACT: The mainobjective ofthis article 
is to present the rates of access to treated 
water in the world and on the American 
continent, as well as to show the importance 
of meteoric water harvesting and encourage 
its use in order to meet the basic needs of 
a large part of the population that does not 
have access. to quality water. For this, the 
method of theoretical documental study was 
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used based on the collection of data made available by large organizations worldwide. The 
results showed agreement on the use of meteoric water, and especially rainwater, for potable 
and non-potable purposes, depending on the use, after physical-chemical treatment. Finally, 
itis concluded that investment in the use of meteoric water can help access to treated water 
for people who do not have access to quality water, as well as alleviate the demand for water 
from public water supply systems, bringing cost savings. for the benefited population. 
KEYWORDS: Rainwater harvesting. Water quality. Water saving potential. 


11 INTRODUÇÃO 


A pandemia da COVID-19 ultrapassou a marca de 600 milhões de casos positivos 
no mundo e já causou quase 6,5 milhões de mortes (WHO, 2022). Ao final do ano de 2019 e 
início de 2020 a Organização Mundial da Saúde (OMS) estipulou medidas preventivas que 
evitam a disseminação rápida da doença e constatou-se que uma das principais medidas 
de contenção do vírus está vinculada à higiene, que é a lavagem das mãos sempre que 
possível. Ou seja, uma das medidas de suma importância para controlar a propagação 
do vírus só poderia ser realizada por aqueles que tivessem o privilégio do acesso à água 
tratada em suas residências e locais de convívio. 

O acesso da população à água potável é um direito básico garantido por lei e esse 
acesso poderia evitar mortes não só por COVID-19, mas também por outras doenças. 
No entanto, de acordo com a UNICEF e a OMS, em 2019 cerca de 3 bilhões de pessoas 
não tinham acesso às instalações básicas de higiene das mãos e cerca de 2,2 bilhões 
nem sequer tinham acesso à água tratada. O saneamento inseguro é um dos maiores 
problemas de saúde pública do mundo. 

A água usada em instalações básicas de higiene deve ser doce e tratada. Essa água 
é captada em mananciais e passa por tratamento para que se torne adequada ao uso da 
população. Porém, existem fontes alternativas de água que podem ser usadas em alguns 
tipos de instalações sanitárias desde que a água seja tratada e se adeque ao uso, como é 
o caso das águas meteóricas. 

As águas atmosféricas são classificadas como águas meteóricas (neblina, névoa 
e precipitação: chuva, neve e granizo) e deposições (orvalho, geada etc.), conforme 
Môller (2008). Dentre as águas meteóricas, a precipitação, sobretudo a água de chuva, 
se destaca por ser de fácil observação, captação e armazenamento. Essa classe de 
água já é aproveitada em diversos países do mundo e aplicada para usos não-potáveis, 
mas demanda certo cuidado em seu manejo e uso. Sendo assim, realizando captação e 
tratamento corretos que ajustem a qualidade da água meteórica para uso em instalações 
sanitárias, é possível ampliar o acesso à água tratada no mundo. 


21 OBJETIVOS 


Este artigo tem como objetivos ressaltar a importância do uso das águas meteóricas 
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e mostrar como o aproveitamento dessa classe de águas pode ampliar o acesso à água 
tratada contribuindo para a adaptação e o enfrentamento do ser humano a um novo 
ambiente imposto pelas mudanças climáticas. 


31 METODOLOGIA 


Foram utilizados os métodos de pesquisa exploratória e explicativa, com base em 
dados recentes disponibilizados em plataformas de grandes organizações mundiais, como: 
Organização Mundial da Saúde (OMS), Organização das Nações Unidas (ONU), Our World 
in Data e Fundo de Emergência Internacional das Nações Unidas para a Infância (UNICEF). 
Essas plataformas auxiliaram no levantamento de dados antes e durante a pandemia, 
como: taxa de acesso à água potável no mundo, índices pluviométricos na américa do 
sul, taxa mundial de mortes associadas ao consumo de água não-potável, taxa mundial de 
pessoas sem instalações seguras de saneamento, região da américa com menor Índice de 
Desenvolvimento Humano (IDH) e acesso a instalações seguras de água tratada. Buscou- 
se, também, conhecer experiências reais existentes no Brasil do uso de água meteóricas, 
com ênfase na água de chuva para o consumo humano, o tipo de tratamento aplicado para 
o seu aproveitamento e seus efeitos econômicos na sociedade local. 


41 RESULTADOS E DISCUSSÃO 


De acordo com o Fundo de Emergência Internacional das Nações Unidas (UNICEF, 
2019) e a Organização Mundial da Saúde (OMS), em 2019 uma a cada três pessoas 
não tinham acesso à água potável do mundo e cerca de 75% da população dos países 
menos desenvolvidos não tinham instalações básicas para higiene e lavagem das mãos. 
Segundo o estudo Global Burden of Disease e a plataforma Our World in Data, as mortes 
relacionadas ao saneamento ocupam 13º lugar em 28 posições e chegaram a alcançar 
mais de 1.2 milhões de mortes no mundo em 2019. A região do continente africano se 
destaca negativamente apresentando os maiores números de mortes (acima de 50 a 100 
por 100 000 pessoas), enquanto os países do continente americano apresentaram entre O 
a 25 mortes por 100 000 pessoas. Em 6 dos 10 países menos desenvolvidos do continente 
americano mais de 40% da população não tem acesso a água tratada, sendo eles Belize, 
Guiana, Guatemala, Nicarágua, Honduras e Haiti. WHO (2020). Dentre os 6 países com 
menores taxas de acesso, a Guiana se destaca negativamente com apenas 27% da 
população com acesso à água tratada, localizada na América do Sul. Em contrapartida, 
apesar de apresentar baixa quantidade de acesso à água de qualidade, a América do Sul 
apresenta altos índices pluviométricos e grande disponibilidade hídrica, como mostra a 
Figura 1. 


Engenharia sanitária e ambiental: Recursos hídricos e tratamento de água 3 Capítulo 3 33 


quútmca qmsos 10N4 

1400 

qq 1300 

' 1200 

1054 1100 

iso RA 1000 

BOLivia Brasito 900 

pts Pas 2051 800 

à Ene) 700 

E rs Camarão s 600 

res f sen. 500 

enc netas mata poco 400 

era Aros mecenaçes 4054 300 

| 200 

À 100 
+ * 5084 
q B0S 

BOM ow 40w 


Figura 1. Mapa geográfico e precipitação acumulada na América do Sul em 2022 (março-setembro). 
Fonte: NOAA (2022). 


De acordo com a Figura 1, os índices de precipitação acumulada são maiores na 
região norte da América do Sul, onde está localizado parte ou inteiramente os países Brasil, 
Peru, Colômbia, Venezuela, Guiana, Suriname e Guiana Francesa. Sendo assim, é possível 
que o baixo acesso a instalações de saneamento seguras para higiene nessa região possa 
ser suprido por fontes alternativas de água tratada, como é o caso da captação das águas 
meteóricas. Tornando, particularmente as águas de chuva, uma importante fonte de água 
para a prevenção da COVID-19 (Kanno et al., 2021). 

De acordo com Pizzo et al. (2021), as águas pluviais devem ser captadas, tratadas 
e armazenadas antes de seu uso. A captação é feita através de um sistema de calhas 
instaladas no telhado de residências e edifícios, tubulações, filtros e reservatórios, 
assim quando ocorre uma precipitação a água escoa por todo o sistema e fica reservada 
aguardando tratamento para ser utilizada, pois apresenta impurezas advindas das telhas 
que teve contato. A qualidade da água é determinada pela análise de alguns parâmetros 
que devem seguir um padrão para cada finalidade de uso, sendo eles principalmente pH, 
cor, turbidez e coliformes. Além disso, fatores naturais também influenciam na qualidade 
da água captada, como tipo, duração e intensidade de chuva, vegetação do local, poluição 
do ar e estação do ano. Para que a água seja usada em fins potáveis, ela deve seguir os 
padrões de potabilidade estabelecido por lei, como é o caso da Portaria GM/MS nº 888/ no 
Brasil (Brasil, 2021). 

Segundo Calheiros et al. (2014), O tratamento de águas pluviais deve prever pelo 
menos a remoção de sólidos, o ajuste do pH e a desinfecção. Geralmente, o tratamento 
é composto por filtração e desinfecção. A filtração tem como objetivo principal a remoção 
de partículas em suspensão na água coletada. A primeira água coletada chamada de first 
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flush é a água que carrega a maior quantidade de impurezas, e por isso comumente é 
descartada. Os filtros são constituídos de camadas de areia e cascalho, e conforme a 
água passa por eles os poluentes ficam retidos. Alguns tipos de material filtrante podem 
ser empregados para complementar o tratamento, como é o caso do carvão ativado que 
consegue agir removendo pequenas partículas orgânicas que influenciam no sabor, odor 
e cor da água. A desinfecção é fundamental para conferir segurança no uso da água. Para 
isso, pode-se utilizar o tratamento por radiação ultravioleta que atua em comprimento de 
onda específico inativando os micro-organismos potencialmente patógenos (Tran et al., 
2021). Outra opção é a desinfecção por cloração e basicamente se resume à adição de 
produtos químicos a base de cloro que agem inativando organismos que podem causar 
riscos à saúde humana. Sendo assim, garantindo que a água passe por todos esses 
processos com os devidos cuidados, é possível que seja aplicada para fins potáveis e não 
potáveis de acordo com a demanda. 

No Brasil, têm-se verificado a implantação de vários sistemas de aproveitamento das 
águas de chuva em áreas isoladas, centros urbanos, indústrias e instituições de ensino, 
ações incentivadas por políticas públicas governamentais para diminuir as desigualdades 
na disponibilidade hídrica e no acesso ao sistema público de abastecimento de água. 
Por exemplo, no município brasileiro de Itapemirim no Espírito Santo, um sistema de 
captação de água de chuva foi implantado em diversos pontos da cidade e apresentou 
retorno positivo para a população, conferindo uma economia de 50% da água cedida pelo 
sistema público de abastecimento de água potável (Bandeira e Stradiotti, 2016). Tornando 
o sistema público mais resiliente frente aos efeitos das mudanças climáticas (Imteaz et al., 
2021). O que foi constatado por Tavares et al. (2022) que há potencial socioeconômico para 
captação de água de chuva em vários municípios brasileiros, sendo influenciado pela tarifa 
e consumo de água. 

A ação de maior destaque é o programa um milhão de cisternas (P1MC) cujo objetivo 
é fornecer água a população rural do semiárido brasileiro (Doss-Gollin et al., 2016). Nesta 
região marcada por fenômenos de seca, Silva et al. (2022) acreditam que o potencial do 
uso de águas de chuva deve ser analisado localmente, pelo tipo de edificação, regime de 
chuvas e custo do serviço público de água. Dentro desta perspectiva, surge no estado 
brasileiro do Pará, o Projeto Amana Katu da Enactus Brasil que disponibiliza sistemas 
de coleta de água de chuva para comunidades quilombolas vulneráveis da Amazônia, e 
também em comunidades rurais de Vazante e Paracatu em Minas Gerais, com o objetivo de 
proporcionar melhores condições de vida e acesso à água potável para a população, visto 
que nessas regiões o acesso à água de qualidade é quase nulo. No sistema desenvolvido, 
a água passa por um filtro que retém sólidos, em seguida passa por um separador que faz o 
descarte da primeira água e a água restante pode ter duas finalidades: a) Fins não-potáveis 
como limpeza; e b) fins potáveis como consumo humano, neste caso a água passa por 
cloração e filtro de carvão ativado, a fim de respeitar aos padrões ambientais (Neiva, 2019). 
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51 CONCLUSÃO 


Por fim, após análise dos dados e exemplos de programas colocados em prática no 
Brasil com o objetivo de levar água potável aqueles que não têm acesso e reduzir a demanda 
de água tratada das concessionárias, conclui-se que a expansão do aproveitamento das 
águas meteórica poderia auxiliar com grande impacto positivo no acesso à água tratada 
em diversas regiões do mundo onde a população não têm acesso à água de qualidade, 
pois ter água disponível não é sinônimo de qualidade. Portanto, medidas que auxiliem na 
redução da taxa de pessoas sem acesso a instalações seguras de água são necessárias, 
visto que este é um dos maiores problemas sanitários do mundo, e captação de águas 
meteóricas, principalmente águas pluviais, tem grande potencial de resolução em parte 
deste desafio. Além de diminuir a pressão pela demanda de água sobre sistemas públicos 
de água proporcionando economia de custo. 
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GESTÃO DAS ÁGUAS TRANSPOSTAS PELO PROJETO SÃO 
FRANCISCO - PISF/ EIXO LESTE: IMPACTOS SOBRE A 
DISPONIBILIDADE HÍDRICA E ESTIMATIVAS DAS TARIFAS 
DE ÁGUA PARA O CANAL ACAUÁ-ARAÇAGI, LOCALIZADO 
NO BAIXO CURSO DO RIO PARAÍBA 
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RESUMO: Tratase de análises dos 
impactos decorrentes do aporte da vazão 
do Eixo Leste do Projeto de Integração 
do rio São Francisco (PISF) sobre as 
disponibilidades do canal Acauã-Araçagi 
e, estimativas das tarifas de água a ser 
cobrada aos usuários deste sistema adutor. 
Principal obra hídrica da Paraíba, o canal, 
concebido em concreto, seção trapezoidal e 
capacidade para vazão máxima de 10 m/s, 
inicia-se no reservatório Acauã, no baixo 
curso do rio Paraíba, estendendo-se por 
112 km até alcançar o leito do rio Araçagi, 
na bacia de mesmo nome. Neste estudo 
foram realizadas simulações utilizando-se 
o modelo ModSIM P32 para avaliações 
das disponibilidades hídricas para o canal, 
nos cenários sem e com a transposição 
das águas do São Francisco. No cenário 
com transposição considerou-se um aporte 
hídrico inicial do projeto de 0,5 mº/s, com 
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acréscimo anual de mesmo valor, até o limite 
da vazão firme de 4,2 mº/s estabelecido 
para o referido eixo. Os impactos sobre os 
valores dos Custos Médios da Água e das 
Tarifas de Água foram analisados para as 
situações demandas para o canal de 0,2 
m/s e de 0,5 m*/s, por meio de cálculos do 
Valor Presente Líquido (VPL) aplicado para 
taxa de retorno de 12% a.a. e horizonte de 
projeto de 30 anos. A variação dos custos 
operacionais e das Tarifas de Água, em 
função das vazões aduzidas, mostrou a 
importância de um planejamento prévio das 
demandas que serão requeridas por parte 
dos estados beneficiados pelo PISF. Estes 
estados devem estar, até a entrada em 
operação do sistema, com suas entidades 
operadoras de suas 
hídricas 


infraestruturas 
instaladas e operantes, com 
seus sistemas de cobrança implantados 
e os valores tarifários definidos, de forma 
a garantir ressarcimento à Operadora 
Federal, os custos operacionais relativos ao 
atendimento das vazões demandadas. 
PALAVRAS-CHAVE: Transposição, outorga 
de uso da água, custos operacionais. 
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ABSTRACT: These are analyses of the impacts resulting from the contribution of the flow 
of the Eastern Axis of the São Francisco river basin and Northern Northeast river basins 
Integration Project (PISF) on the availability of the Acauã-Araçagi channel and estimates of 
the water tariffs to be charged to users of this adductor system. The channel is designed in 
concrete, trapezoidal section and capacity for maximum flow of 10 m3/s. It starts in the Acauã 
reservoir, on the low course of the Paraíba River, extending for 112 km to the Araçagi River. 
In this paper, simulations were performed using the ModSIM P32 model to evaluate the water 
availability for the canal, in the scenarios without and with of the water of the São Francisco 
river. In the scenario with transposition, an initial water contribution of the project of 0.5 m/s 
was considered, with an annual increase of the same value, up to the limit of the firm flow 
rate of 4.2 mº/s established for the said axis. The impacts on the values of the Average Water 
Costs and Water Tariffs were analyzed for the situations demands for the channel of 0.2 m%/s 
and 0.5 m/s, through calculations of the Valor Presente Líquido (LPV) applied for a retum 
rate of 12% a.a. and project horizon of 30 years.The variation in operating costs and water 
tariffs, due to the flows inducted, showed the importance of prior planning of the demands that 
will be required by the states benefited by the PISF. These states must be, until the system 
has come into operation, with their operating entities of their installed and operating water 
infrastructures, with their collection systems deployed and the tariff values defined, in order to 
guarantee reimbursement to the Federal Operator, the operational costs related to meeting 
the required flows. 

KEYWORDS: Transposition, granting water use, operating costs. 


INTRODUÇÃO 


O Projeto de Integração de Águas do Rio São Francisco com Bacias do Nordeste 
Setentrional (PISF) prevê uma entrada de água, via leito natural do rio Paraíba até atingir 
o açude Poções no município de Monteiro. A adução terá uma parcela de 4,2 mº/s de 
vazão firme e 5,8 m*/s de vazão flutuante a ser disponibilizada, dependendo das condições 
hídricas do reservatório de Sobradinho. Dentre as regiões que serão beneficiadas pelo 
referido eixo, situam-se: o Cariri, o Curimataú e o Seridó, nas quais é observado o maior 
problema concernente ao abastecimento de água do Estado da Paraíba, e onde está 
inserido o município com menor índice pluviométrico do Brasil: Cabaceira, com uma média 
anual em torno de 300 mm. 

A escassez de água é agravada pelo déficit crônico do balanço hídrico regional, que 
envolve a pluviometria e a evaporação. O clima tropical, com elevada temperatura e alta 
insolação, ocasiona taxa de evaporação de até 2.500 mm ao longo do ano, estabelecendo 
sempre a uma situação de escassez hídrica na região. O manejo das reservas hídricas 
é refém da incerteza quanto à afluência hídrica para o próximo período de chuva. Dessa 
maneira, a operação dos reservatórios é limitada em função da necessidade de se 
armazenar água como garantia para anos vindouros, quando poderá ou não haver uma 
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recarga hídrica. 

No contexto de distribuição das águas transpostas pelo Eixo Leste do PISF no rio 
Paraíba, foi projetado o canal Acauã-Araçagi, partindo do reservatório Acauã, no baixo 
Curso do rio Paraíba, estendendo-se por 112 km estendendo-se até alcançar o reservatório 
Araçagi, na bacia do rio Mamanguape. Com capacidade para transportar até 10 m3/s tem 
por finalidades: abastecimento humano de 33 municípios situados na Zona da Mata e no 
Agreste paraibano, com população estimada em 590 mil habitantes e demanda de 0,938 
mê/s estimada para 2035), dessedentação animal e irrigação de uma área potencial de 16 
mil hectares. 

Com objetivos de avaliar a sustentabilidade hídrica e de estimar as respectivas 
tarifas de água a serem cobradas aos usuários do canal adutor, com vista ao ressarcimento 
dos custos à Operadora Federal, foram realizadas simulações do sistema hídrico do rio 
Paraíba para os cenários sem e com incremento das vazões aduzidas pelo Eixo Leste/ 
PISF, considerando duas situação de demandas hídricas a serem atendidas pelo referido 
canal. 


O EIXO LESTE DO PISF 


O Eixo Leste do Projeto de Integração do rio são Francisco com bacias do Nordeste 
Setentrional (PISF) inicia-se com a captação no reservatório de Itaparica, em Floresta, no 
estado de Pernambuco até desaguar no açude Poções, em Monteiro, no estado da Paraíba. 
A interligação do Eixo Leste do PISF com a Paraíba se dará através de uma derivação do 
reservatório Barro Branco (trecho V), passando pelo portal de entrega de água (PBO1L), 
na divisa dos estados da Paraíba e Pernambuco, até desaguar no açude Poções, na bacia 
hidrográfica do rio Paraíba (Lima, 2011). Com extensão aproximada de 290 km, desnível 
geométrico de 304 m e capacidade para transportar a vazão máxima de 28 mº/s, tem como 
objetivo principal atendimento às demandas hídricas das bacias do rio Paraíba e dos rios 
Moxotó e Ipojuca, no estado de Pernambuco (FUNCATE, 2001). O eixo é composto por: 
13 barragens, 211 km de segmentos de canal em concreto, seção trapezoidal, 2,5 km de 
aquedutos e 22,5 km de túnel (figura 1). 
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Figura 1 — Integração do Eixo Leste/PISF com o rio Paraíba 


O SISTEMA HÍDRICO ESTUDADO 


A Bacia do Rio Paraíba 


A bacia hidrográfica do Rio Paraíba, localizada ente as latitudes 6º51'31” e 8º26'21” 
Sul e as longitudes 34º48'35” e 37º2'15” a Oeste de Greenwich é a segunda maior do 
Estado da Paraíba, abrangendo uma área de 20.071,82 Km?, que representa 38% do 
território e uma população de 1.828.178 habitantes, correspondente a 52% da população 
total. Considerada uma das bacias mais importantes do semi-árido nordestino, ela é 
composta pela sub-bacia: Taperoá, Alto Paraíba, Médio Paraíba e Baixo Paraíba, a jusante 
do reservatório Acauã até sua foz no oceano Atlântico. 

Em termos de climatologia, de acordo com a classificação de Kôeppen, a região 
do Alto Curso do Rio Paraíba possui clima do tipo BSwh' (semi-árido quente) e do Baixo 
Curso do Rio Paraíba vigora o clima do tipo Aw' (úmido). A vegetação natural predominante 
é caatinga hiperxerófila, hipoxerófila, floresta caducifólia e subcaducifólia na região do Alto 
e Médio cursos. Na região do Baixo Curso existem algumas áreas com a vegetação nativa 
da Mata Atlântica e ecossistemas associados. As precipitações médias anuais variam de 
400 mm, na região do Alto Paraíba até 1700 mm na região do Baixo Paraíba. A evaporação 
média anual varia entre 1.300 mm e 2.500 mm. Os solos mais representativos são: Bruno 
não Cálcico, Litólico, Solonetz Solodizado, Regossolo e Cambissolo. 
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Canal adutor Acauã-Araçagi 


O canal adutor Acauã-Araçagi é uma obra em concreto, seção trapezoidal, extensão 
aproximada de 112 Km, com capacidade para transportar a vazão de 10 m*/s, para operação 
por gravidade. Com vista ao aproveitamento dos recursos hídricos advindos do Eixo Leste 
do PISF, a obra foi projetada com as finalidades: abastecimento humano de 33 municípios 
situados na Zona da Mata e no Agreste paraibano, com população estimada em 590 mil 
habitantes e demanda de 0,938 mº/s (horizonte 2035), dessedentação animal e irrigação 
de uma área potencial de 16 mil hectares. Principal obra hídrica do estado da Paraíba, 
orçada aproximadamente em 934 milhões de reais, inicia-se no açude Acauã, no médio 
curso do rio Paraíba, estendendo-se até alcançar a barragem de Araçagi, na bacia do rio 
Mamanguape, no qual está previsto desaguar 2,5 mº/s. A outorga da vazão para o canal foi 
estabelecida pela Agência Executiva de Gestão das Água do Estado da Paraíba (AESA) 
em 2,0 m?/s. A figura 2 mostra o sistema hídrico estudado: alto, médio e baixo cursos do rio 
Paraíba e o canal adutor Acauã-Araçagi. 


q 
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druvdo, $ 
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y , 
at ms Regiões do Alto, Médio e Baixo Curso do 
Rio Paraiba 
Figura 2 — Sistema hídrico estudo: Bacia do rio Paraíba/Canal Acauã-Araçagi 
METODOLOGIA 


Para avaliação das disponibilidade hídricas para o canal adutor foram realizadas 
simulações do sistema hídrico do rio Paraíba utilizando-se o modelo ModSIM P32 (Labadie, 
1989), para os cenários sem e com aporte hídrico advindo do Eixo Leste do PISF para o 
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canal, por meio de entrega da vazão firme de 4,2 mº/s no reservatório Poções. A modelagem 
deste sistema hídrico está mostrado na Figura 3. 


“4 


Figura 3 — Esquema para simulação do sistema hídrico do rio Paraíba 


As estimativas dos custos da água decorrentes dos aportes hídricos para o canal 
adutor tiveram como referência os custos operacionais relativos à formulação do sistema 
tarifário do PISF (FGV, 2005), base para as emissões da Outorga de uso da água (ANA, 
2005a) e do CERTOH do empreendimento (ANA, 2005b). As estimativas e suas implicações 
sobre as tarifas de água, a serem cobradas dos usuários do canal adutor, contemplaram os 
três conjuntos de valores dos custos operacionais: custos fixos, custos variáveis e BDI da 
Operadora Federal (Lima, 2013). As mesmas foram obtidos por meio do cálculo do Valor 
Presente Líquido (VPL). 


VPL=SJ [(REC-CO)/(1+i)] (1) 


Em que: 
VPL — Valor Presente Líquido (US$); 
REC — Receitas auferidas da cobrança pelo uso da água (US$); 
CO — Custos operacionais envolvidos nos sistemas hídricos (US$); 
i— Taxa de retorno (%); e 
n- Horizonte de projeto (anos) 
Os Custos Médios da Água (CMA) foram definidos como a relação ente os custos 
totais operacionais do sistema hídrico estudado (PISF, sistema hídrico do rio Paraíba e 
canal Acauã/Araçagi) e as vazões disponibilizadas ao sistema hídrico estudado, mostrado 
na Figura 2. A Tarifa de Água (TA) foi definida como o menor valor para obtenção do VPL 
positivo, considerando o horizonte de projeto igual a 30 anos e taxa de retorno de 12% a.a. 
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Cenário sem a transposição (C1) 


No cenário sem aporte hídrico, as simulações do sistema hídrico do rio Paraíba 
foram realizadas priorizando as demandas de abastecimento dos reservatórios integrantes 
do mesmo, inclusive, a vazão ecológica à jusante do reservatório Acauã, estabelecida em 
0,5 mº/s. As vazões máximas disponibilizadas para o canal adutor foram obtidas por meio 
de simulações com a série de 56 anos de vazão pseudo-histórica. Não foram considerados 
aportes hídricos de outras sistemas para o canal. Para as estimativas das tarifas de água 
somente os custos de operação e de manutenção do canal foram analisadas, considerando 
uma taxa de incremento anual de 5%. O custo relativo à gestão a Operadora Estadual não 
foi considerado, pois a mesma não havia sido definida oficialmente. 


Cenário com a transposição (C2) 


No cenário com a transposição, as simulações do sistema hídrico estudado foram 
realizadas com periodicidade anual, considerando aportes hídricos do PISF, iniciando com 
a vazão de 0,5 mº*/s e incremento anual de mesmo valor, até atingir a vazão firme a ser 
fornecida pelo Eixo Leste para a bacia do rio Paraíba, igual a 4,2 m/s. Os valores das 
tarifas a serem praticados, para o ressarcimento dos custos à operadoras Federal, foram 
obtidos por meio do cálculo do Valor Presente Líquido (VPL). Sobre os custos operacionais 
do canal adutor foram adicionados: os custos de operação e manutenção, gestão, outorga 
e energia, relativos as disponibilizações das vazões requeridas. 

Para cada vazão disponibilizada pelo PISF ao sistema hídrico estudado foram 
obtidas as vazões máximas disponibilizadas pelo reservatório Acauã para o atendimento as 
demandas do canal Acauã-Araçagi. Nas simulações foram consideradas as perdas hídricas 
em trânsito, entre os reservatórios Poções (entrada da vazão transposta) e Acauã (tomada 
d'água do canal adutor), para obtenção das respectivas vazões a serem bombeadas pelo 
projeto. 

Os custos relativos a demanda e ao consumo de energia foram obtidos considerando 
o bombeamento de 21 horas, em dias úteis, e de 24 horas nos finais de semana e feriados, 
utilizando-se as tarifas “fora de ponta”. Os valores praticados neste estudo foram: 3.25 
US$/Kw e 73.63 US$/Mwnh, para as tarifas de demanda e consumo, respectivamente. 

Os custos relativos ao processo de disponibilização das vazões do Eixo Leste para o 
sistema estudado foram repassados para o cálculo das tarifas de água para o canal adutor, 
proporcionalmente ao ganho de vazão, conforme apresentado na Equação 2. 

CEC = (GCQ 


/QP, )xCTE (2) 


disp disp 
Em que: 
CEC — Custo relativo ao consumo de energia para o canal; 


GCQ,..— Ganho de vazão para o canal, come semo aporte hídrico datransposição; 


disp 
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QP 
CTE — Custo total de energia para bombeamento da Q 


asp — Vazão disponibilizada pelo PISF para o sistema hídrico do rio Paraíba; e 


me dO PISF para o sistema 
PB. 

Neste estudo, o custo relativo à outorga foi obtido proporcionalmente à vazão firme 
máxima a ser disponibilizada para o estado da Paraíba, cujo valor definido para o Eixo 
Leste/PISF foi de US$ 1,107.87 por ano. Para ambos os cenário, também foram realizadas 
simulações considerando a situação de aumento da demanda de 0,2 mº/s (S1) para 0,5 
mê/s (S2), com incrementos anuais de mesmos valores. Foram realizadas análises dos 
impactos sobre os custos médios da água e sobre a tarifa de água a ser praticada ao 


usuários do canal adutor para as duas situações de demanda. 


ANÁLISES DOS RESULTADOS 


No cenário sem atransposição (C1), a vazão máxima disponibilizada pelo reservatório 
Acauã para o canal adutor foi de 0,95 mº/s, que corresponde à vazão regularizada com 
100% de garantia. No Quadro 1 estão os valores obtidos dos Custos Médios da Água 
(CMA), das Tarifas de Água (TA) e os horizonte de atendimento das demandas pelo canal, 
nas duas situações estudadas. Na situação de demanda S1, o horizonte de atendimento foi 
de 2 anos e a tarifa de água foi inferior em 0.007 US$/mº (18,5%), em relação a situação 
S2, cujo horizonte obtido foi de 5 anos. Este resultado mostra que, neste cenário, quanto 
mais rápida for utilizada a vazão disponibilizada pelo reservatório Acauã para o canal 
adutor, menor será o valor da tarifa de água. Este fato é decorrente das maiores receitas 
auferidas com a cobrança da água, devido a vazão ser consumida no menor horizonte de 


atendimento a demanda. 


Demanda Vazão Valores (US$/m?) Horizonte 
(canal) (m/s) CMA TA (Ano) 
S1 0,20 0.0385 0.0376 5 
S2 0,50 0.0385 0.0306 2 


Quadro 1. Custo Média da Água e Tarifas de Água para as situações de demandas S1 e S2 


No Quadro 2 estão mostrados os custos operacionais envolvidos na disponibilização 
dos aporte hídricos do PISF ao sistema estudado, envolvidos nas análise dos impactos 
sobre os valores dos custos médios da água e da tarifa de água, para as duas situações de 
demanda no cenário com a transposição (C2). 
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Custos operacionais médios anuais 
St (mi US$) S1(%)  S2(miUS$) S2(%) 


Item dos custos 


Operação e manutenção (canal Acauã-Araçagi) 1.153 15,7% 1.153 9,6% 
Operação e manutenção (reservatórios-PB) 0.045 0,6% 0.064 0,5% 
Operação e manutenção (infraestrutura-PISF) 1.967 26,7% 2.802 23,2% 
Gestão (PISF) 0.047 0,6% 0.080 0,7% 
Demanda de energia 0.218 3,0% 0.369 3,1% 
BDI da Operadora Federal 0.335 4,5% 0.487 4,0% 
Outorga de uso da água 0.206 2,8% 0.350 2,9% 
Consumo de energia 3.032 41,2% 6.028 50,0% 
Taxa da Operadora Federal 0.364 4,9% 0.723 6,0% 
Total 7.367 100,0% 12,057 100,0% 


Quadro 2. Custos operacionais anuais da infraestrutura do sistema hídrico para as situações de 
demandas S1 e S2. 


Os dados mostram que os maiores custos operacionais da infraestrutura dos 
sistemas hídricos envolvidos no estudo estão relacionados ao consumo de energia, com 
valores anuais de 3.032 milhões de dólares e de 6.028 milhões de dólares, correspondente 
aos percentuais 41,2% e de 50% dos custos totais, para as situações de demandas S1 e 
S2, respectivamente. 

Para as vazões disponibilizadas pelo PISF, variando de 0,5 mº/s, com incremento 
anual de mesmo valor até a vazão firme de 4,2 m/s, foram obtidos os valores dos Custos 
Médios da Água (CMA) e das Tarifas de Água (TA) correspondentes a situação de demanda 
do canal S1 (0,20 mº/s). Os valores estão apresentados no Quadro 8. 

Os custos médios da água diminuíram e as tarifas de água aumentaram com os 
acréscimos das vazões disponibilizadas pelo PISF. Até a vazão disponibilizada pelo projeto 
de 1,5 mº/s, os valores das tarifas de água foram inferiores aos custos médios da água. 
Para a vazão disponibilizada igual a 2,0 mº/s, o custo médio da água e da tarifa de água 
apresentaram a menor diferença de valores, igual a US$ 0.003. Também se observa 
que, para a referida vazão, os valores CMA e TA, bem como a relação entre custo de 
bombeamento e os custos operacionais, foram os que mais se aproximaram dos valores 
médios da série de vazão disponibilizada pelo projeto. Este fato pode ser indicativo de que 
a partir desta vazão, os custos relacionados ao consumo de energia (Quadro 2) passaram 
a ser preponderantes sobre os demais custos operacionais do sistema hídrico estudado. 
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Vazão 


agent Crhemu  Joume Cuso cua turio (bem 
(m/s) 
PISF Canal (Canal-PISF) Canal (Hm?) ia a (%) 
050 1,31 0,36 o 4418 2826 0196 0064 224 
100 1,45 0,50 0,14 28830 3779 0131 0078 319 
150 1,61 0,66 0,16 43254 497 114 004 390 
200 1,78 0,83 0,17 57672 6201 ot os! 44,1 
250 1,92 0,97 0,14 72090 7395 0103 0125 482 
300 201 1,06 0,09 86508 8396 0097 0140 518 
350 2,14 1,19 0,13 100926 9624 0095 0156 545 
420 242 1,47 0,28 21 11861 0098 18! 56,8 
Média 0.118 0.119 43,6 


Quadro 3. Vazões disponibilizadas pelo PISF e respectivos CMA e TA para situação (S1) 


Neste cenário foram também analisados, para cada aporte hídrico do PISF, os 


aumentos dos valores das tarifas de água e dos custos relativos aos ganhos de vazão 


disponibilizada pelo canal adutor, na situação de demanda S1. Os resultados no quadro 


4 mostram que, para o incremento de 1 mº/s na vazão disponibilizada do canal adutor 


custará aproximadamente 7.920 mi US$ por ano. Para o bombeamento da vazão firme 


estabelecida para o Eixo Leste (4,2 m*/s), que proporcionará ao referido canal adutor um 


ganho de vazão de 1,47 mº/s, custará anualmente o valor de 11.861 mi US$, não incluídos 


os custos relativos aos serviços da Operadora Estadual. Observa-se ainda que, para a 


vazão disponibilizada pelo PISF de 2,5 mº/s, a relação custos operacionais totais e ganho 


de vazão disponibilizada apresentou o valor mais próximo da média da série de vazão, 


igual a 7.623 mi US$ por ano. 


EE 
(m/s) (m/s) (mi US$m'/s) (miUS$ano/m'/s) (US&m) (USSm) (%) | (US$m?) 
0,5 0,36 2.826 7.850 0.064 -— == 0.026 
1,0 0,50 3.779 7.559 0.078 0.015 22,8 0.041 
1,5 0,66 4.927 7.465 0.094 0.016 20,6 0.057 
2,0 0,83 6.201 7.471 0.111 0.017 17,8 0.074 
2,5 0,97 7.395 7.623 0.125 0.014 12,4 0.087 
3,0 1,06 8.396 7.921 0.140 0.015 12,3 0.103 
3,5 1,19 9.624 8.087 0.156 0.015 10,8 0.118 
4,2 1,47 11.861 8.068 0.181 0.025 16,1 0.143 
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Média 6.876 7.755 0.119 0.017 16,1 0.081 


Quadro 4. Valores das tarifas e dos custos unitários da vazão disponibilizada pelo canal 


O menor aumento do valor da tarifa, em temos absolutos, foi observado quando se 
aumentou a vazão disponibilizada de 2,0 m*/s para 2,5 m?/s, com o valor igual a 0.014 US$/ 
mê. Entretanto, em termos percentuais, a menor variação foi observada entre as vazões 
disponibilizadas: 3,0 mº*/s e 3,5 mº/s, com valor igual a 10,8%. Em termos de valores 
médios, correspondentes ao aumento da tarifa de água para o canal, observa-se que o 
valor absoluto obtido para a vazão disponibilizada de 2,0 mº/s se aproxima da média e, o 
valor médio percentual de aumento da tarifa (16,1%) foi observado para a vazão 4,2 m/s, 
que corresponde a máxima disponibilizada pelo Eixo Leste do PISF. 

No Quadro 5 estão apresentados os valores obtidos para os Custo Média da Água 
(CMA), para as Tarifas de Água (TA), os respectivos horizontes de atendimento para as 
duas situações de demandas do canal Acauã-Araçagi, no cenário com a transposição (C2). 


Demanda Vazão Valores (US$/mº) Horizonte 
(canal) (ms) CMA TA (Ano) 
S1 0,20 0.0979 0.1806 18 
S2 0,50 0.0995 0.1152 5 


Quadro 5. Custo Média da Água e Tarifas de Água para as situações de demandas S1 e S2 


Os resultados mostram que o aumento da vazão do canal em 0,30 m*/s, na situação 
de demanda S2, causou um acréscimo de 1,6% no valor dos Custos Médios da Água, e 
um decréscimo de US$ 0.0655 (36%) no valor da Tarifa de Água, em relação a situação 
de demanda S1. Outra observação importante, refere-se a diferença significativa entre os 
valores da TA e dos CMA, na medida em que maiores vazões foram disponibilizadas pelo 
PISF para o canal adutor. Este fato é decorrente das próprias definições destes parâmetros. 
Enquanto o valor dos CMA está associado aos custos operacionais e aos volumes a serem 
disponibilizados pelo projeto, a TA está relacionada aos custos operacionais do projeto 
e das infraestruturas hídricas internas dos estados beneficiados pelo projeto, no caso, o 
Estado da Paraíba. Neste cenário, assim como observado no cenário sem a transposição 
(C1), os resultados mostraram que menores valores da Tarifa de Água (TA) foram obtidos 
para menor horizonte de atendimento as demandas do canal Acauã-Araçagi. No sistema 
estudado, redução do horizonte de 13 anos para 5 anos. 


CONSIDERAÇÕES FINAIS 


No cenário sem transposição (C1), as disponibilidades hídricas do canal, a partir 
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da vazão regularizada do açude Acauã, não teve grandes impactos sobre os valores dos 
CMA e das TA, os quais se mantiveram muito próximos. Este fato está associado aos 
menores horizontes de atendimento as demandas do canal adutor e, principalmente, por 
não incluírem os custos de energia para o bombeamento das vazões do PISF. No cenário 
com aporte hídrico do PISF (C2), os resultados mostram que, para o incremento de 1 mº/s 
na vazão disponibilizada pelo canal adutor custará aproximadamente 7.920 mi US$ por 
ano. A vazão firme estabelecida para o Eixo Leste (4,2 mº/s), que proporcionará ao canal 
adutor um ganho de vazão de 1,47 mº/s, custará anualmente o valor de 11.861 mi US$, não 
incluídos os custos relativos à gestão da infraestrutura hídrica pela Operadora Estadual. 
Neste cenário, os valores dos CMA e da TA foram próximos na média da série de vazões 
disponibilizadas pelo Eixo Leste do PISF. Nos dois cenários estudados, os resultado 
mostraram que quanto mais rápida for utilizada a vazão disponibilizada pelo projeto ao 
canal adutor, menor será o valor da Tarifa de Água. Este fato é decorrente das maiores 
receitas auferidas com a cobrança da água, devido a vazão ser consumida no menor 
horizonte de atendimento à demanda. A variação dos valores das Tarifas de Água, em 
função das vazões aduzidas do PISF, mostrou a importância de um planejamento prévio 
das demandas que serão requeridas por parte dos estados beneficiados pelo projeto. Neste 
estudo, ressalta-se a importância de que, até a entrada em operação do canal adutor, 
o Estado da Paraíba deve estar com suas entidades operadoras de suas infraestruturas 
hídricas instaladas e operantes, com seus sistemas de cobrança implantados e com os 
valores tarifários definidos, de forma a garantir o ressarcimento à Operadora Federal, os 
custos operacionais relativos ao atendimento das vazões demandadas do projeto. 
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RESUMO: O arsênio é um elemento tóxico, 
sendo o ecossistema aquático vulnerável 
as fontes de contaminação antrópicas e de 
biodisponibilidade natural. O monitoramento 
de arsênio em baixas concentraçõesaliado 
aos procedimentos de remediação de 
águas se faz cada vez mais necessário para 
promover um controle eficiente junto às 
normatizações ambientais. Neste trabalho 
foi proposto o uso da argila expandida 
como um adsorvente para remover 
arsênio do meio aquoso, aplicado em um 
dispositivo de extração miniaturizado em 
ponteira descartável. A extração ocorre de 
forma reversível com a dessorção, sendo 
obtidas condições ótimas para detecção do 
arsênio em sistema de geração de hidretos. 
A proposta pode ser considerada uma 
alternativa viável, que permite recuperar 
e enriquecer a quantidade de arsênio 
extraído, abrindo novas possibilidades 
para ser usado como reagente de segunda 
geração. 

PALAVRAS-CHAVE: Arsênio, remediação, 
argila, adsorvente. 
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USE OF A MINIATURIZED DEVICE FOR REMOVING ARSENIC IN AQUEOUS 
ENVIRONMENTAL USING EXPANDED CLAY AS EXTRACTION ADSORBENT 
ABSTRACT: Arsenic is a toxic element, and the aquatic ecosystem is vulnerable to human 
sources of contamination and natural bioavailability. The monitoring of arsenic in low 
concentrations combined with water remediation procedures is increasingly necessary to 
promote an efficient control along with environmental regulations.In this work it was proposed 
the use of expanded clay as an adsorbent to remove arsenic from the aqueous medium, 
applied in a miniaturized extraction device in a disposable tip.The extraction occurs reversibly 
with desorption, obtaining optimal conditions for the detection of arsenic in a hydride generation 
system.The proposal can be considered a viable alternative, which allows recovering and 
enriching the amount of extracted arsenic, opening new possibilities to be used as a second- 

generation reagent. 
KEYWORDS: Arsenic, remediation, clay, adsorbent. 


11 INTRODUÇÃO 


O arsênio é um elemento reconhecido por apresentar alta toxicidade mesmo presente 
em baixas concentrações em amostras provindas de matrizes ambientais e alimentícias. 
Essas amostras, uma vez contaminadas podem integrar as vias de absorção no organismo 
animal, e acumular espécies tóxicas de arsênio nos tecidos fisiológicos, comprometendo a 
regulação metabólica (MIHUCZ, 2022). 

O arsênio se manifesta na natureza em suas formas inorgânicas, principalmente 
nas espécies trivalente, como arsenito (As**) e pentavalente, como arsenato (As**).Em 
menor proporção o arsênio se manifesta em formas orgânicas, como espécies metiladas 
(MMA e DMA), e associadas em açúcares ederivados metabólitos oriundos de vitaminas 
e amoninoácidos, tais como arsenocolina (AsC) e arsenobetaína (AsB).Esses últimos são 
denominados compostos organoarsenicais, e sáomais abundantes no ecossistema marinho. 
Compostos orgânicos de arsênio também são utilizados como aditivos em defensivos 
agrícolas, com por exemplo o Roxarsone? e Daconate(BJORKLUND et al., 2020). 

As formas inorgânicas de arsênio são cerca de 100 vezes mais tóxicas em relação 
aos compostos orgânicos, sendo que a espécie trivalente é considerada a mais crítica, 
sendo 60 vezes mais tóxica em comparação com a sua forma mais oxidada (CHATTERJEE 
et al., 1995). Compostos organoarsenicais, não apresentam toxicidade significativa para a 
biota aquática. 

O ambiente aquático é mais eficiente em acumular e biodisponibilizarespécies de 
arsênio do que no solo (LITER; MORGADA, 2009). Recursos hídricos minerais e águas 
oceânicas próximas a regiões de depósitos auríferos, de ferro, manganês, e que sofrem 
influência de atividade vulcânica são mais vulneráveis a contaminação por arsênio. Essa 
contaminação pode acarretar na acumulação nociva de arsênio em alguns alimentos 
cultivados sob a forma de irrigação, como o arroz (PAULELLI et al., 2019; SANTOS et al., 
2021). 
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No Brasil, é relatado um histórico de contaminação de arsênio na região do 
Quadrilátero Ferrífero em Minas Gerais, devido a atividades de exploração de jazidas 
de ouro, dos quais depósitos de pirita (FeS), são comumente encontradas (BORBA; 
FIGUEIREDO; CAVALCANTI, 2004). Relatos críticos decontaminação emáguas 
influenciadas pela atividade de mineração, também foram reportados na cidade mineira de 
Paracatu (REZENDE et al., 2015). 

A exposição humana a elevados níveis de arsênio ocorre principalmente pelo 
consumo de água contaminada e por alimentos preparados com essas águas. A Agência 
de Proteção Ambiental Americana (EPA-US), em conformidade com a Organização Mundial 
da Saúde (OMS), revisou o limite tolerado de arsênio de 50 para 10 mg L”, em finalidades 
de classificação da potabilidade de água (WHO, 2010). O CONAMA também atendeu a 
essa nova deliberação, reforçando os valores orientadores de qualidade da água através 
da Resolução 357, de 17 de março de 2005, e pela Portaria 2914, de 12 de dezembro de 
2011 do Ministério da Saúde. 

Diante dessa problemática, torna-se necessário o investimento em tecnologias de 
monitoramento do teor de arsênio em águas naturais, e quando constatado o grau de 
contaminação, disseminar procedimentos que visem a remediação ou mitigação dos níveis 
tóxicos das espécies de arsênio disponíveis no meio aquático. Procedimentos de redução 
da concentração de arsênio através da sua remoção por adsorção tem recebido grande 
destaque, devido as vantagens de simplicidade de operação, eficiência e baixo custo. Além 
disso, pode ser enquadrado como uma alternativa ecologicamente viável, quando se usam 
adsorventes provindos do reuso adequado de materiais naturais, como os de constituição 
lignocelulósica (PANIAGUA; COSTA; COELHO, 2021) 

A adsorção de arsênio foi satisfatoriamente investigada através do uso da farinha da 
casca de banana (PANIAGUA et al., 2021) e da casca de maracujá (SANTOS et al., 2021), 
visando aplicação no tratamento de águas contaminadas. Tratando-se de materiais de 
constituição mineral, o uso de argilas apresenta características físico-químicas compatíveis 
para serem agregadas às vantagens da classe de adsorventes naturais, possibilitando seu 
estudado nos procedimentos de remoção de arsênio em água. 

A grande maioria dosprocedimentos experimentais baseados em adsorção, são 
realizados através de testes em batelada. Se tratando de resíduos aquosos contaminados, e 
da expectativa de concentração de arsênio a níveis traço,ainda são poucos os trabalhos que 
podem ser efetivamente aplicados em uma escala menor.Nessas condições, normalmente 
se empregam pequenos volumes de amostra, logo o uso de dispositivos miniaturizados 
para procedimentos de adsorção se torna mais adequada, uma vez que podem ser menos 
dispendiosas e menos onerosas do que as são executadas em batelada. A extração em 
ponteira descartável conhecida como procedimento DPX é uma alternativa que tem se 
destacado ultimamente para o desenvolvimento de procedimentos de adsorção em fase 
sólida, que pode ser adequada em escala miniaturizada (BORDIN et al., 2016; PINTO; 
QUEIROZ, 2015). 

No procedimento DPX, uma pequena quantidade da fase sólida adsorvente é 
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acondicionada no interior de uma ponteira de pipeta descartável, do qual é acoplada em 
uma seringa para promover uma rápida dispersão de ar, no momento em que se aspira 
a amostra. A extração do analito é prontamente realizada pela dispersão eficiente da 
amostra nas partículas do adsorvente, promovendo um contato dinâmico e instantâneo. O 
procedimento DPX pode ser facilmente adaptado a condições de infraestrutura laboratorial 
disponível, e compatível para determinações moleculares e elementares (FERNANDES et 
al., 2014; SILVA et al., 2022; TOMASIN et al., 2021). 

Neste trabalho, é proposto o uso da argila expandida como umafase sólida adsorvente 
a ser empregada em um dispositivo de extração de arsênio através dos mecanismos de 
adsorção, utilizando uma adaptação do procedimento DPX, com o objetivo de ser aplicado 
a remoção de arsênio em amostras aquosas que podem estar contaminadas. 


21 PROCEDIMENTO EXPERIMENTAL 


2.1 Materiais e Reagentes 


Soluções de trabalho de As** na concentração de 20,0 ug L'! foram preparadas 
por diluição adequada da solução estoque (1000 mg L!) em água deionizada. A solução 
estoque foi obtida pela dissolução de NaAsO, em meio de HCI 0,3 mol L”. 

Para o sistema de extração, bem como na etapa de análise e determinação de 
arsênio, foram preparadas soluções aquosas de HCI partindodo reagente concentrado de 
grau P.A. Soluções do redutor NaBH foram estabilizadas em meio de NaOH 0,5% (m/v). 

Soluções de NaOH e HCl 0,1 mol/L foram utilizadas para ajustar o pH do meio 
aquoso das amostras. O ajuste foi feito através de um pHmetro de bancada com eletrodo 
de membrana de vidro combinado com eletrodo de referência Ag/AgCI. 


2.2 Dispositivo de extração em ponteira descartável 


Para montagem do sistema DPX adaptado para uma escala miniaturizada, foram 
utilizadas ponteiras plásticas descartáveisde micropipetas automáticas de 5 a 10 mL (tipo 
Gilson). Para promover a aspiração e a dispersão da amostra, foram utilizadas seringas 
plásticas com a abertura injetável ajustada ao encaixe com a base superior da ponteira. Na 
extremidade inferior da ponteira, foi introduzida uma pequena porção de lã de vidro, afim 
de promover filtragem e limpeza durante a aspiração da amostra. 

A argila expandida, utilizada como adsorvente extrator foi obtida em casas de 
jardinagem comerciais. Partes maiores do sólido foram primeiramente quebradas com 
auxílio de um martelo e em seguida trituradas em moinho e peneiradas. As partículas de 
argila trituradas foram lavadas com água deionizada e secas em estufa. Foram pesadas 
40,0 mg do pó da argila expandida e introduzidas manualmente para dentro da ponteira. 
Em seguida, foi feita uma lavagem com 1,0 mL de água deionizada para condicionamento 
do sólido e posterior inserção do filtro superior com lã de vidro para apoio no encaixe com 
a seringa. 

As etapas de extração de arsênio consistiram na aspiração de uma alíquota de 4,0 
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mL da solução amostra ou solução teste de As!'20,0 mg L', seguida da aspiração de ar 
para promover a dispersão com o adsorvente durante 40 s. Após esse tempo de contato, a 
solução aspirada foi descartada. O conjunto dessas etapas foi denominada como um ciclo 
de extração, do qual ocorre-se a remoção de arsênio da matriz aquosa. 

A fim de se obter a recuperação da quantidade de arsênio adsorvida na argila, 
realizou-se uma etapa de eluição. Essa etapa consistiu na aspiração de uma alíquota de 
0,5 mL da solução eluente de HCl para dessorver o arsênio da fase sólida, seguida da 
aspiração de ar (aproximadamente 3,0 mL). Posteriormente, a solução foi dispensada em 
frasco de leitura da amostra para análise. 

As etapas do procedimento foram ilustradas na Figura 1, representando o ciclo de 
extração e eluição, executados no dispositivo de DPX adaptado: 


<> 
H 


Condicionamento Aspiração Aspiração Dispensar 


Eluição 
amostra dear Amostra 


Figura 1: Ilustração das etapas experimentais da extração de arsênio por DPX. 
Fonte:: Adaptado de BORDIN et al., 2016 


2.3 Sistema analítico de determinação de arsênio e otimização 


O arsênio foi detectado e determinado nas amostras aquosas por Espectrometria de 
Absorção Atômica por Chama com Geração de Hidretos (HG-AAS). Um sistema comercial 
de Geração de Hidretos VGA 77 (Varianº, Victoria, Austrália), foi acoplado ao espectrômetro, 
equipado com atomizador de cela de quartzo posicionada sobre o queimador e alinhado ao 
feixe óptico emitido da lâmpada de catodo oco de arsênio, operando nas especificações do 
fabricante. A representação esquemática do sistema VGA 77, com desing de funcionamento 


Engenharia sanitária e ambiental: Recursos hídricos e tratamento de água 3 Capítulo 5 


55 


em sistema de fluxo contínuo, está representada pela Figura 2. 


Tubo atomizador de quartzo 


Bomba Peristáltica Bobina de Reação 


Amostra 


HCl /H20 


Separador gás-líquido 


No (100 mL min? 


Figura 2. Representação esquemática do sistema comercial de geração de hidretos VGA77. 
Fonte: Manual do usuário (VGA 77) 


A absorbância foi tomada como sinal analítico, da qual foi otimizada em função 
da recuperação de arsênio obtida após eluição da quantidade adsorvida na fase sólida 
do dispositivo. A otimização foi realizada de modo multivariado, através da construção da 
matriz de planejamento fatorial 2º, das quais foram avaliadas as seguintes variáveis: |) 
concentração de HCl; ii) concentração do redutor NaBH,; iii) vazão da amostra; iv) vazão 
do redutor NaBH, e v) vazão do HCI. Para um refinamento da triagem realizada pelo 
planejamento fatorial, um modelo de superfície de resposta foi desenvolvido por meio do 
Planejamento Composto com Ponto Central (CCD). O tratamento estatístico da otimização 
dos dados analíticos foi executado com auxílio do software Statistica 7.0 (StatSoft). 

Os níveis estudados para cada variável referente aos planejamentos fatorial e CCD 
estão representados na Tabela 1. 


Variáveis 
Nível HCl NaBH, Vazão amostra Vazão HCI Vazão NaBH, 
(mol L-?) Yo(mAv) (mL min!) (mL min?) (mL min?) 


Planejamento Fatorial 


-1 1,0 0,10 0,8 0,8 0,8 

0 5,5 0,45 1,1 2,3 2,3 

+ 10,0 0,80 1,4 4,5 4,5 
Planejamento Composto Central (CCD) 

-1,68 5,0 0,1 go Es 1,5 
41 6,0 0,5 -— -— 2,3 
0 7,0 0,8 -— -— 4,5 
+1 8,0 1,2 = = 5,5 
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+1,68 9,0 1,5 -—- = 6,7 


Tabela 1. Níveis das variáveis investigadas na otimização multivariada para construção das matrizes 
dos planejamentos fatorial 2º e CCD. 


31 RESULTADOS E DISCUSSÕES 


Inicialmente foi investigada o efeito da variação do pH do meio aquoso da amostra 
sobre a adsorção de As** na argila expandida. Os valores de pH testados, variaram na 
faixa de 2,0-10,0. Nesse ensaio, a solução amostra que foi aspirada, deixada em contato 
com o adsorvente e posteriormente dispensada pelo dispositivo DPX (sobrenadante),foi 
levada para análise no sistema HG-AAS. A absorbância foi usada como sinal analítico 
para determinação qualitativado teor de arsênio remanescente. A remoção foi calculada 
mediantea diferença relativa da absorbância da solução inicial (As!! 20,0 ug L'), pela 
absorbância do sobrenadante. Os resultados são apresentados na Figura 3A. 

Uma vez avaliada a influência do pH sobre a extração, foi realizado um estudo 
da capacidade de recuperação do arsênio removido. O ensaio foi baseadona dessorção 
do arsênio provindo das mesmas soluções empregadasno ciclo de extração,em relação 
aosseusrespectivos valores de pH. O As* foi dessorvido da argila expandida por eluição 
com solução de HCI 0,1 mol L'!. A recuperação foi expressa a partir do ganho de sinal 
obtido pela etapa de eluição, a fim de que o procedimento seja promissor para efeito de 
pré-concentração em amostras contaminadas a níveis traço de arsênio. O ganho de sinal 
foi calculado pela razão entre as absorbâncias da solução eluída em relação ao sinal da 
solução de As!! 20,0 ug L!, tomada como controle. A Figura 3B apresenta os resultados 
obtidos da recuperação de arsênio. 


(A) (B) 
100.00 7.00 
97.78 6.65 
95.00 94.15 Rue 
Ra — 5.00 
se q 
o 90.00 9113 > 400 0.35 
o 90.52 E 
(6) “a 
É 85.00 q 3:00 0.05 
e E 
e 200 
80.00 8810 E 0.50 
O 1.00 
0.28 
75.00 0.00 
O 2 4 6 8 10 12 O 2 4 6 8 10 12 
-1.00 
pH pH 


Figura 3. (A) Efeito da variação do pH sobre a remoção de Asº* em meio aquoso por adsorção na argila 
expandida. (B) Estudo da influência do pH sobre a recuperação de Asº** mediante etapa de eluição 
(Controle: 0,0496 de absorbância para As** 20,0 ug L"!). 


Fonte: Os autores (2022) 
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Os resultados indicam que a argila expandida apresenta boa capacidade de adsorção 
de As** do meio aquoso, com remoção média de 92,34%. Esse valor representa um 
desempenho satisfatório das condições de extração empregadas no sistema DPX. Melhor 
resultado de remoção pôde ser alcançado em pH 4,0. Nesse pH, a superfície da argila se 
encontra carregada positivamente, uma vez que a partir desse estudo, foi verificado um 


valor de pH de carga zero próximo a neutralidade (pH, ..= 7,2). Nessa condição de acidez, 


PCZz 
o As!!, se manifesta na forma do ácido arsenioso (H,AsO,; pKa, = 9,2; pKa, = 12,1), como 
uma espécie eletricamente neutra, uma vez que se encontra completamente protonada. 

O favorecimento da remoção de arsênio nessas condições, dão indícios de que 
o mecanismo de adsorção pela argila não ocorre de forma específica, por interações 
químicasnem eletrostáticas. A retenção pode ser razoavelmente explicada através da 
formação da membrana de Donnan, como acontece nos processos de troca iônica por 
exclusão. Em amostras aquosas, os grupos funcionais da argila se tornam hidratados, sendo 
permeada por uma camada negativamente carregada (membrana de Donnan),que permite 
que somente moléculas não dissociadas possam atravessá-la (EITH et al., 2006).Nesse o 
caso, como o H.AsO, é uma molécula neutra, esta consegue atravessar a membrana de 
Donnan e ser adsorvida pela argila através do estabelecimento de ligações de hidrogênio. 

A eficiência do ciclo de extração em pH 4,0 se deve ao fato de que a espessura 
da camada de hidratação delimitada pela membrana de Donnan se torna maior em 
comparação com o que ocorre em pH 2,0 e 6,0. Esse fenômeno permite que um ligeiro 
excesso de H,AsO, seja adsorvido pela argila. 


O aumento da adsorção observado em pH 10,0 (>pH,.,; >pKa,), pode ser justificado 


PCcz” 
em função da contra-carga (positiva) da atmosfera iônica que circunda a espécie aniônica 
H,AsO,.Devido ao aumento da força iônica do meio provocada pela adição de NaOH, pode- 
se promover interações de natureza eletrostática favoráveis na superfície do adsorvente, 
que nessa condição encontra-se com sítios ativos carregados negativamente. 

Já com relação a recuperação de arsênio, os resultados indicam de forma geral, 
que a solução de HCI na concentração de 0,1 mol L não foi suficiente para promover 
uma eluição significativa. Observa-se um maior ganho de sinal para arsênio recuperado na 
solução amostra com pH 2,0. Sugere-se que nessa condição fortemente ácida, a superfície 
da argila pode estar saturada em sítios de carga positiva, ocorrendo o pareamento iônico 
direto com ânions cloreto. Dessa forma, a membrana de Donnan pode não ser eficientemente 
estabelecida durante a extração, o que facilita a dessorção de arsênio. 

Diante do que foi exposto, no intuito de melhorar a eficiência da etapa de eluição e 
recuperação do arsênio, foi realizado um estudo variando a concentração da solução de 
HCl na faixa de 0,10-1,25 mol L'!. Para esse estudo, foi fixado o valor de pH 4,0 durante o 
ciclo de extração. Os resultados são mostrados na Figura 4 a seguir. 
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Concentração (mol L'?) 


Figura 4. Estudo da variação da concentração da solução de HCI como eluente para recuperação de 
arsênio por dessorção (Obs.: Agua deionizada foi tomada como amostra em branco. O teste sem DPX 
representa a leitura da solução controle de As** 20,0 ug L'?). 


Fonte: Os autores (2022). 


A concentração de HCI 0,25 mol L! foi escolhida como a melhor condição para 
aprimorar a recuperação do arsênio durante a eluição. Essa concentraçãorepresenta 
um ganho de sinal de 48,33 vezes em relação ao sinal obtido da solução controle, que 
não passou pelo procedimento de extração DPX.Concentrações maiores do eluente HCl 
passam a contribuir mais para a degradação da rede cristalina mineral da sílica, reduzindo a 
dimensão das partículas de argila dispersas, e aumentando consequentemente a superfície 
de contato do adsorvente. Isso faz com que a dessorção do arsênio seja dificultada, 
resultando numa queda na capacidade de recuperação do arsênio. 

Amelhora do sinal alcançado para o arsênio recuperado na etapa de eluição se torna 
primordial para promover boa sensibilidade de detecção no sistema HG-AAS. Visando a 
aplicação em amostras de águas naturais, espera-se que os prováveis níveis de arsênio 
oriundos por fontes de contaminação se enquadrem em uma ordem de grandeza muito 
baixa. Para tanto, é fundamental certificar que a técnica analítica escolhida para detecção, 
apresente capacidade para quantificar concentrações a níveis por partes por bilhão (ppb, 
equivalente a ug L!) com precisão e exatidão aceitáveis. 

A sensibilidade de detecção de As!!! pelo sistema HG-AAS foi checada através da 
construção de uma curva de calibração (Figura 5), do qual pôde-se verificar que os limites 
de detecção e quantificação obtidos pelo equipamento foram respectivamente 2,52 e 8,42 
Hg L respectivamente. 
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q 
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Concentração As!" (ppb) 


Figura 5. Curva de calibração obtida para padrões aquosos de Asº* no sistema HG-AAS. 


Fonte: Os autores (2022). 


Afim de garantir o comprometimento da detectabilidade de As** a níveis traço, frente 
a possibilidade de pré-concentração no momento em que se proporciona a recuperação do 
arsênio extraído pela argila, foram feitos ensaios de otimização multivariada das variáveis 
que influenciam o sistema HG-AAS. Os efeitos de cada variável estudada, bem como os 
efeitos de interação conjunta (de segunda ordem) entre as variáveis estão apresentadas no 
gráfico de Pareto (Figura 6). Os resultados foram tratados ao nível de 95% de confiança, 
dos quais podem ser considerados significativas as variáveis que apresentaram valor de 
efeito que ultrapassaram a linha de corte em p=0,05. 

De acordo com o Gráfico de Pareto, a concentração do redutor NaBH, foi a variável 
que apresentou maior efeito significativo, sendo a de maior importância no sistema HG- 
AAS, seguido da concentração do HCI. Os efeitos positivos dessas variáveis indicam que 
maiores concentrações do redutor e do ácido favorecem o sinal analítico. A reação de 
geração do hidreto volátil de arsênio (AsH,) é promovida pela ação do redutor em meio 
suficientemente ácido, para que haja a evolução de hidrogênio nascente, reativo com a 
espécie trivalente de arsênio. 
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(2) [NaBH4] + 

Curvatr. 

(1) [HCl] 

2by4 

4by5 

1by2 

1by5 

1by3 

3by4 

2by3 

(3) Vazão Amostra | 0,6789026 

3by5 -0,382756 

2by5 + -0,375305 

1by4 -0,334329 

(5) Vazão HCl 0,3212913 
(4)Vazão NaBH4 0,1164099 


p=,05 


Standardized Effect Estimate (Absolute Value) | 
Figura 6. Gráfico de Pareto dos efeitos das variáveis obtidas pelo planejamento fatorial 2º. 


Fonte: Os autores (2022). 


Uma das principais vantagens da otimização multivariada em relação ao modo 
univariado, é a possibilidade de investigar a interação entre as variáveis, além do número 
reduzido de experimentos (NETO; SCARMINIO; BRUNS, 1995). Tal fato pode ser 
observado pela importância da interação entre os efeitos das vazões do ácido e redutor, 
que isoladamente não apresentem qualquer efeito significativo sobre o sinal. O efeito 
positivo dessa interação (4'5) indica um comportamento sinérgico, com favorecimento da 
absorbância, quando as soluções do ácido e do redutor são introduzidas no sistema de 
fluxo do HG-AAS nos respectivos níveis maiores de vazão simultaneamente. A vazão da 
amostra não apresentou influência sobre o sinal de detecção, bem como sua interação com 
outras variáveis também não resultou em efeitos significativos. 

Um outro aspecto importante que deve ser levado em consideração, é a contribuição 
da curvatura sobre o modelo matemático linear fornecido pelo planejamento fatorial. O 
experimento realizado nas condições do ponto central (nível 0), é importante para estimar 
o efeito da concavidade sobre o plano do domínio experimental. O gráfico de Pareto indica 
que o efeito da curvatura foi significativo, e que um modelo linear não se mostra adequado 
para uma estimativa das condições ótimas de detecção. Para tanto, um modelo matemático 
quadrático foi desenvolvido a partir da construção de uma matriz CCD (Central Composite 
Desing). 

Para a construção do CCD foram tomadas as variáveis de maior importância 
indicadas pelo gráfico de Pareto, como a concentração do HCl, bem como a concentração 
e vazão do NaBH,.Os novos níveis estudados para essas variáveis no Planejamento 
Composto Central partiram da indicação de triagem revelados pelo planejamento fatorial. 
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A vazão do HCl foi fixada em 4,5 mL min” e não foi incluída na matriz CCD, uma vez que 
experimentos realizados em vazões maiores poderiam ocasionar a formação em excesso 
do gás H, que podem provocar diluição do hidreto volátil de cela atomização de quartzo, 
comprometendo a sensibilidade de detecção, além do risco de ignição. 

O tratamento dos resultados obtidos pela execução da matriz CCD, possibilitou a 
construção das superfícies de respostas quadráticas, dos quais são ilustradas nas Figuras 
7A, Be 7C a seguir: 


Figura 7. Superfícies de resposta obtidas por CCD. (A) [HCI/[NaBH ]; (B) [HCl]'vazão NaBH,; (C) 
[NaBH ]'vazão NaBH,. 


Fonte: Os autores 


A Análise de Variância (ANOVA) indicou a razoabilidade do modelo quadrática, com 
uma variação explicada de 66,27% dos dados, com uma média de erro puro de 6,0'10ºº. 
As superfícies acusaram um ponto crítico que pode ser considerada representativa das 
condições ótimas de detecção pelo sistema HG-AAS, sendo elas: HCl 6,6 mol L! e NaBH4 
0,9% (m/v), operando em uma vazão de 4,7 mL min”. 


41 CONCLUSÃO 


A argila expandida se mostrou eficiente para adsorção de As** em meio aquoso, com 
remoção superior a 90%. O dispositivo de extração por DPX se mostrou uma alternativa 
viável para miniaturização, possibilitando a recuperação de arsênio e seu enriquecimento 
durante a eluição com HCl. Condições ótimas de detecção pelo sistema HG-AAS puderam 
ser estimadas por meio de ferramentas de otimização multivariada. A proposta deste 
trabalho se mostra promissora para remediar recursos hídricos contaminados com arsênio, 
além de abrir novas possibilidades para ser recuperado e usado como reagente de segunda 
geração. 
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